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ABSTRACT: 

This report contains the material presented in the form of panel discus- 
sions at the Thirteenth Annual Air Force Office of Scientific Research Con- 
tractors' Conference on Physical Energetics held at the Naval Postgraduate 
School, Monterey, California, from the 5th to the 7th of May 1970. 

It was the purpose of this conference to facilitate the exchange of 
latest information on research in non-chemical energy sources, release 
mechanisms and conversion processes. The conference was attended by 
the AFCSR contractors and grantees performing research in these areas to- 
gether with representatives of the Army, Navy, Air Force, DCD, NASA and 
AEG. 



COMPILATION OF RESEARCH SUMMARIES 



THIRTEENTH AFOSR CONTRACTORS' MEETING 

ON 

PHYSICAL ENERGETICS 



5-7 May 1970 
Naval Postgraduate School 
Monterey, California 



2 



THIRTEENTH AFOSR CONTRACTORS' MEETING ON PHYSICAL ENERGETICS 

5-7 May 1970 

AGENDA 



Tuesday 
5 May 1970 

0800 Registration 

0830 Opening Remarks NPS and AFOSR Representatives 



0900 


Injection and Plasma 


Panel 1 


(Part 1) 


Page 5 




Generation 








1030 


COFFEE 








1045 


Injection and Plasma 


Panel 1 


(Part 2) 






Generation 








1230 


LUNCH 








1330 


Turbulence and Waves 


Panel 2 


(Part 1) 


Page 54 


1500 


COFFEE 








1515 


Turbulence and Waves 


Panel 2 


(Part 2) 




1700 


ADJOURN 








Wednesday 










6 May 1970 










0830 


Confinement and Instabilities 


Panel 3 


(Part 1) 


Page 72 


1015 


COFFEE 








1030 


Confinement and Instabilities 


Panel 3 


(Part 2) 




1230 


LUNCH 








1330 


Calculations and Instrumen- 


Panel 4 


(Part 1) 


Page 121 




tation 








1500 


COFFEE 








1515 


Calculations and Instrumen- 


Panel 4 


(Part 2) 






tation 








1700 


ADJOURN 








1845 


RECEPTION COCKTAILS 








1930 


BANQUET Professor E. C. 


Haderlie, Guest Speaker 


Thursday 










7 May 1970 










0830 


Low Temperature Plasmas 


Panel 5 


(Part 1) 


Page 145 


1015 


COFFEE 








1030 


Low Temperature Plasmas 


Panel 5 


(Part 2) 




1230 


LUNCH 








1330 


Round Table Discussion of 


Panel 6 


(Part 1) 





Future Plasma Applications 



3 



Thursday (Cont'd) 
7 May 1970 



1500 


COFFEE 




1515 


Round Table Discussion of 


Panel 6 (Part 2) 


1700 


Closing Remarks 


NPS and AFOSR Representatives 


1715 


ADJOURN 





4 



THIRTEENTH AFOSR CONTRACTORS' MEETING ON PHYSICAL ENERGETICS 

5-7 May 1970 

PANEL 1 . INJECTION AND PLASMA GENERATION 



Name 


Topic 


Dr. Willard H. Bennett 
North Carolina State University 
Department of Physics 
Raleigh, North Carolina 27607 


Production of Extremely Dense Hot 
Plasmas 


Professors. Gruber 
Department of Engineering 
Case Western Reserve University 
Cleveland, Ohio 44106 


Controlled Beam -Plasma Interaction 
in Non-Uniform Magnetized Plasmas 


Dr. A. Haught 

United Aircraft Research Labs 
United Aircraft Corporation 
East Hartford, Connecticut 06108 


Laser Irradiated Single Particle 
Plasmas 

V 


Professor C. D. Hendricks 
Department of Electrical Engineering 


Electron Sheath Extraction of Ions 
from a Plasma 


Columbia University 
New York, N Y. 10027 


Thermonuclear Reactions Initiated in 
Particle -Particle Collisions 


Dr. M. Lubin 

Dept, of Mechanics & Aerospace 
Science 

University of Rochester 
Rochester, N. Y. 14627 


Interaction of Laser Produced Plasmas 
with a Magnetic Field 


Professor F. Schwirzke 
Department of Physics 
Naval Postgraduate School 
Monterey, California 93940 


Laser Created Plasma and Collisionless 
Shock Wave Phenomena 


Dr. JohnG. Siambis 
Department of Electrical Engineering 
Carnegie-Mellon University 
Pittsburgh, Pennsylvania 15213 


Streaming Plasma Interaction with 
Longitudinal Magnetic Fields 



5 



PRODUCTION OF EXTREMELY DENSE HOT PLASMAS 

Willard H. Bennett, North Carolina State University, Raleigh, N. C. 

The immediate objective in this program is to develop methods for 
concentrating the beams from very high current high voltage pulsed electron 
beam machines into streams with very small diameters, by using the self- 
magnetic field of a linear pinched discharge. Machines delivering 30 
nanosecond pulses of 100,000 amperes at 5,000,000 volts have been used 
in this work. The method which has been used is to produce a linear 
pinched discharge in argon at 0.2 Torr and while the pinched discharge 
is still of the order of one centimeter in diameter, to project the pulse 
of high energy electrons through a 0. 002-inch metal sheet into the end of 
the linear pinched plasma column so that the self-magnetic field of the 
linear pinch serves to draw the relativistic electron beam towards the 
axis . 

Experimental confirmation has been attained for several unusual 
characteristics predicted theoretically for these relativistic beams as 
they interact with a pinched discharge and which are not familiar with 
ordinary discharges. The injection of the pulsed beams into the pinched 
column induces reverse current densities initially which are everywhere 
equal and opposite to the injected relativistic current densities. It has 
also been predicted that the production of induced reversed current makes 
the residual electrons in the beam channel move backwards so that they 
are unpinched, that is, they are magnetically driven out of the beam channel, 
whereupon the injected relativistic electron beam is able to pinch itself 
far more intensely so that the radius of the relativistic beam is reduced 
to a very much smaller value than it had at injection. It is anticipated 
that by properly adjusting the parameters and by developing methods for 
delivering the relativistic electron beam through the thin metal foil within 
a half-angle of divergence of less than one degree, a resulting radius of 
beam less than one-tenth millimeter should be attained. 

The production of the electron beam within the small angle of diver- 
gence is to be called **narrow injection.'* Methods for producing narrow 
injection in the existing high energy pulsed electron beam machines is 
the princi p al objective of this program at this time and it is anticipated 
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that to attain this it will be necessary to develop some rather sophisticated 
methods for making the pinched discharge interact appropriately with the 
discharge system in the high voltage tube of the high energy machines « These 
methods are under development. 

An ancillary objective of this program to run concurrently with the 
above is to study the processes in targets which will be produced when 
electron beams with power densities of the order of 1,000c 000,000 jOOO, 000 
watts per square centimeter strike various kinds of solid state targets. 

In the early phases of this work we have already found that our basic 
.athode which is slightly concave and approximately two inches in diameter 
can be made to deliver 100,000 amperes at the full 5,000,000 volts across 
a two- inch gap with no damage whatsoever to either the cathode or the anode. 
This demonstrates that the concept of a plasma cathode for these purposes 
is sound and that the physical processes in these discharges can be in- 
vestigated effectively. These are being studied both experimentally and 
theoretically preparatory to proceeding with the investigation of the 
interaction of the pinched discharge with the discharge system in the 
high voltage tube. 

Although the scope of the immediate program is limited to the pro- 
duction of extremely high power densities in beams, there are some very 
important areas of application which will become accessible as soon as the 
objective of the program is attained. A few of these can be mentioned as 
follows s (1) the production of thermonuclear power by nuclear fusion at 
plasma densities more than 10,000 times as great as any ever attained 
before on earth; (2) the production of transuranic elements in quantities 
greater than micrograms rather than atom-by-atom as presently done in the 
biggest of accelerators and hence to produce such elements at very much 
higher atomic numbers than presently attainable; (3) production and study 
during very short intervals of very small bits of plasmas similar to those 
in the interior of stars and in nuclear explosions; and (4) the coherent 
acceleration of small clumps of ions of the order of 100,000,000,000 ions 
to energies greater than 10,000,000,000 electron volts per nucleon (if the 
ions are, for example uranium, energies greater than 1,000,000,000,000 
electron-volts per ion) and projecting these accelerated ion clumps of 
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selected composition into solid state targets of selected composition, to 
produce nuclear processes which never occurred before in the history of the 
universe except possibly at the moment of the '*big bang” if indeed such 
was the origin of the universe. 
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CONTROLLED BEAM-PUSMA INTERACTION IN NON-UNIFORM MAGNETIZED PLASMAS 

by 

. '*■ ■- 

S. Gruber 

Division of Electrical Sciences and Applied Physics 
Case Western Reserve University 

Oscillations in beam-plasma systems are usually generated from 
uncontrolled fluctuations in space charge. The observed spectra of such 
oscillations are dominated by the frequencies exhibiting maximum growth 
rates and are generally broad band as the result of multimode excitation 
from these fluctuations. Recently efforts were made to control beam-plasma 

n 21 

instabilities by modulation of the electron beam'- ’ ^ and preparation of 
its distribution function. These techniques have been sufficiently 
successful to give impetus to more experimental efforts using beam-plasma 
interactions in a controlled manner. We have undertaken an experimental 
program to study controlled beam-plasma interactions in an inhomogeneous 
medium with these aims : 

(1) "Plugging” Mirror Machines 

(2) Acceleration of Plasmas 

A large and successful effort is being directed toward the rf accelera- 

rsi 

tion of plasma at the Centre d' Etudes Nucleaires de Saclay in France. '■ ■' 

Their efforts are concentrated on experiments where rf fields are imparted to 
a plasma (which may be also produced by these same fields) in a gradient of 
magnetic field so that the electron cyclotron frequency matches the wave 
frequency in the gradient. This results in the acceleration of plasma into 
the low field region. The rf energy is either applied in a resonant struc- 
ture or in a travelling wave structure. The experimental program to follow 
is designed to investigate the effect of localized beam-plasma interactions 
in a gradient of static magnetic field upon the acceleration and heating of 
plasma. 

Experiment 

The electron beam is density modulated by an rf source at the frequency 
where the wave interaction is to be studied. This is accomplished by the 
experimental set up shown in Fig. 1. The beam cathode (on the left) is 
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mounted flush witli the interior wall of a re-entrant cavity. The grid is 
mounted flush witli the other inside wall and is at the same time insulated 
from the wall to allow for a grid bias. The grid is approximately 0.3 mm 
from the cathode. Tliis small spacing is necessary in order that transit time 
effects will be negligible up to 500 MIlz. The anode is about 7.5 mm in front 
of the grid. The cavity is insulated from ground and is run at a negative 
potential with respect to the grounded anode. 

The emerging modulated beam is injected through a region of corkscrew 
magnetic field, which is immersed in a region of axially directed uniform 
magnetic field. The corkscrew effects a non- adiabatic conversion of part of 
the axial velocity of the beam electrons into rotational motion. The beam 
enters into the mirror region where it is to interact with the plasma present. 
A second and very much weaker beam is injected from the other side of the 
mirror. It is axially directed and unmodulated and is used as a method of 
probing the fields within the confining region. A collector for the probing 
beam is located to the left of the mirror. 

The electron density as well as the magnetic field varies along the 
axis of the macliine. If the energy of the beam electrons on emerging from 
the corkscrew field is properly chosen the beam could be nearly reflected from 
the mirror maximum. This would have the effect of enhancing the growth rate 
of tlie instability. In Fig. 2 we show a sketch of the plasma density, np , 
axial magnetic field strength, and beam density, n^^ , as a function of 

the distance from the electron gun toward the mirror maximum. 
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LASER IRRADIATED SINGIE PARTICIE PLASMAS 



A, F. Haughty D. H. Polk^ W. J. Fader ^ and J. C. Woo 

United Aircraft Research Laboratories 
East Hartford, Connecticut 

The production in vacuum of isolated, high density, high temperature plas- 
mas and the investigation of the confinement of these plasmas in simple mirror 
and minimum-B magnetic fields are being studied under AFOSR support and in 
parallel programs supported by Corporate and USAEC funds. This technique has 
been demonstrated to be an effective means of producing a spherically-expanding, 
uncontaminated plasma at high density and temperature.^ Using the low fre- 
quencies available with a recently developed, variable -frequency (lO-UoO Hz) 
suspension voltage source, particles up to radius can be conveniently 

suswnded to yield plasmas containing up to 10^"^ atoms and, with the existing 
ruby laser pulse (7 joule, I 5 nsec), having an average energy of 200 eV . 

Mass spectrometer measurements show that the plasma contains an indetectable 
fraction (<2 percent) of impurities and, at plasma energies in excess of 100 eV, 
consists of H"^ and predominantly triply ionized lithium ions. 

These plasmas can be produced within established magnetic containment 
fields, and calculations examining the plasma formation and expansion in the 
presence of a magnetic field show that the plasma heating and initial expansion 
are not significantly affected for plasmas formed within magnetic fields up to 
tens of kilogauss field strength. As a result the plasma energy is the same as 
that obtained for a free expansion plasma, and the plasma expansion is initially 
spherically symmetric. However, after the heating phase of the plasma develop- 
ment, interaction of the plasma with the magnetic field results in stopping of 
the cross-field plasma motion and thermalizat ion of the expansion kinetic 
energy to produce a high temperature trapped plasma. Experimental studies con- 
firm the results of these magnetohydrodynamic calculations and, in addition, 
provide information on the subsequent containment of the trapped plasma by 
simple mirror and minimum-B magnetic fields. Charge collection measurements 
show that the cross-field plasma expansion is suppressed by a magnetic contain- 
ment field and that more than 99 percent of the plasma leaves the trapping region 
along magnetic field lines, i.e., within the magnetic field mirror loss cones. 
Comb probe measurements permit determination of the transverse dimensions of 
the contained plasma and show that, throughout the plasma decay from the magne- 
tic field, the plasma radius does not increase appreciably beyond that calculated 
for the initial trapping under the conditions of the experiment. This result 
implies that, on the time scale of the present experiments, instabilities do not 
act to produce significant cross-field movement of the plasma. Expansion 
velocity measurements along the direction of the loss cone show the increase in 
axial velocity expected for an expansion which has been altered from a spherical 
to a one -dimensional (axial) form. This result implies both that the transverse 
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expansion has been stopped and that the trapped plasma is thermalized upon 
capture by the magnetic field. Further evidence for plasma thermalization^ 
predicted theoretically as consequence of shock collisional interactions in the 
magnetic field trapping process^ is given by line radiation studies showing 

reduced recombination and by time resolved retarding potential probe measurements 
of the ion energy distribution late in the plasma density decay. Retarding 
potential probe measurements show a contained ion temperature, equal to - 50 
percent of the laser plasma average energy over the 50-200 eV plasma energy 
range examined. Microwave measurements of the plasma density decay show that, 
in sufficiently large minimum-B magnetic containment fields, the plasma decay 
time is equal to that predicted for Coulomb collisional scattering of a hydrogen 
plasma into the magnetic field mirror loss cones. Further support for a Coulomb 
scattering loss mechanism of the plasma density decay is provided by the 
dependence of the plasma decay time over the factor of four temperature range 
examined. Plasma decay times shorter than predicted by Coulomb scattering are 
observed in simple mirror fields and with lower strength minimum-B fields; in 
each case these are correlated experimentally with a large random cross -field 
plasma loss suggestive of a flute instability which is apparently suppressed for 
sufficiently large minimum-B field gradient. 

In the present research period, time dependent mass spectrometer studies 
of the decaying plasma have been carried out and show rapid escape of the highly 
ionized lithium ions in the plasma followed by a more gradual decay of the 
hydrogen in agreement with a Coulomb scattering loss mechanism. These measure- 
ments also establish that the contained plasma late in its decay is composed 
primarily of hydrogen, a result used in comparison of the measured plasma decay 
time with that predicted on the basis of Coulomb scattering. Flux coil measure- 
ments have been performed which show that the spherical plasma expansion is 
stopped and the plasma trapped at /3 1, as predicted by a simple one -dimensional 

model of the plasma magnetic field interaction. Initial plasma capture at ^ 1 

appears inconsistent with subsequent Coulomb decay of a mirror confined plasma 
and has led to a study of the limitations on laser plasma capture in open-ended 
confinement field geometries. From the results of this study, the plasma capture 
and containment are separate processes, and their different dependence on tem- 
perature leads to a maximum plasma temperature for optimum long-term containment .5 
At low temperatures, the field diffuses rapidly into the plasma, but the life- 
time is limited by the large Coulomb scattering loss. Increasing the plasma 
temperature to obtain a longer Coulomb scattering time also increases the 
characteristic time for the magnetic field to diffuse into the plasma. From 
the dependence of these two times on plasma temperature, theoretical calculations 
show that there is an optimum temperature of 1 keV which results in the greatest 
amount of plasma captured for the longest time . 

However, both the temperature and the subsequent Coulomb decay time of a 
1 keV plasma are too small for significant fusion output, and laser heating 
cannot be used to directly produce a magnetically confined fusion plasma. This 
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result^ however^ does not restrict the utilization of a laser produced plasma 
as a seed plasma from which a thermonuclear plasma is then generated within the 
magnetic confining field. In this approach^ the laser produced plasma is used 
as a target for charge exchange with more energetic neutrals which are then 
trapped by the confining magnetic field. The high captured plasma density 
attainable with a laser produced plasma makes this a particularly attractive 
scheme^ and the use of a neutral injection beam permits replacement of particles 
which have been lost by decay from the plasma. Calculations of the charge ex- 
change heating and beam ionization have been carried out and show that a lOO mA^ 
20 keV neutral injection beam incident upon a 1 keV laser produced plasma of 
10^^ particles can result in the production of a 20 keV plasma at a density of 
2x 10^3 cm"3^ a plasma of definite thermonuclear interest. A most attractive 
feature of this approach is that lasers to produce a 1 keV target plasma and 
neutral injection systems to provide the required 100 mA beam current lie within 
the pre sent state -of -the -art . 

On the basis of these results^ a neutral beam injection-laser produced 
target plasma program has been undertaken. An existing 100 G¥ peak power ^ 

0.1 nsec pulse duration laser will be used to produce the initial target 
plasma. Ion beam sources^ which were developed under AFOSR support in another 
program (Contract AF33(65T)-12385 have been obtained and will be adapted to 
provide neutral atom beam for injection into the target plasma. The AFOSR 
supported studies will be concerned with the neutral beam injection plasma 
heating which represents a logical continuation of the current program and is a 
prerequisite to the attainment of a steady state laser produced thermonuclear 
plasma. Studies of the confinement and stability of the plasma will be carried 
out in a program supported in part by the U.S. Atomic Energy Commission. 
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r:LECTRON SHEATH EXTRACTION OF IONS FROM A PLASMA 



C. D. Hendric s. A. VV. Dipert and S. Grodzinsky 



Charg('d Particle Research Laboratory 
Department of Electrical Engineering 
University of Illinois 



Introduction 

It has been proposed to create a source of ions by passing a sheet 
beam of electrons over a plasma o The sheet beam will act essentially as a 
negative electrode and will thus extract positively charged ions from the 
plasma while repelling the electrons from the plasma o The resulting beam 
of ions may then be focussed and used as an ion source to be used in fields 
such as electric propulsion and controlled thermonuclear research. 

There are several requirements that an ion source should satisfy to 
be useful in electric propulsion » In order to obtain the highest thrust 
possible, it is necessary that the ion beam be of a high current density. 

Also the ion beam should have ideally only one species of ion in order that 
the efficiency in the use of the propellant be made as large as possible. 

To have a uniform thrust, it is obvious that the beam must also have a 
uniform density. Ion sources made for other uses will, of course, have 
other requirements dictated by their particular uses and needs. 

Several different kinds of ion sources have been developed. Among them 
are the contact cesium-tungsten source, the duo-pla smatron , and the magnetron 
source. In addition, W. D. Owen(l) spoke of extracting ions from a plasma 
by using an electrode at a negative potential. The contact cesium-tungsten 
source was utilized by R. X, Meyer (2,3,4) who used a sheath of electrons 
to extract the ions. Magnetron sources have also been examined by several 
authors. Redhead (5) used this device for various pressure measurements 
while Jepsen and Helmer (6) studied various characteristics of both the 
normal and inverted magnetrons. Masek (7), in his investigations of the 
inverted magnetron as a bombardment ion engine, discovered that the radial 
ion density in this configuration was much more uniform than in the normal 
magnetron configuration (by a factor of 5). This would, of course, lead to 
a more uniform ion beam upon extraction from this source. Cobic and Matic 
(8) examined several characteristics of gas ionization in a magnetron and 
noted the effect of the magnetic field upon ionization. In addition, Cobic 
et al . (9) used a magnetron to produce a new kind of ion source for high- 
melting point materials. Gabovich (10) described several different gas 
discharge sources, including the duopla smatron and plasmatron, and also 
studied such basic problems as the ionization of plasmas and the extraction 
of ions from plasmas. 

Space charge limitations as expressed by the Langmuir-Childs law are 
the most severe limitations in most ion sources. The law states that: 



where J is the current density, V is the applied voltage, d is the anode 
cathode separation and k is the quality factor. There are several ways the 
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the current density in this relation may be increasedo The separation d, 
if decreased with the voltage V held constant, will obviously increase J» 

A major problem which exists when a solid electrode is used for extraction 
is that if the separation distance, d, is made too small, disruptive vol- 
tage breakdown may occur. This is one of the major reasons that a sheet 
electron beam is to be used for extraction instead of a solid electrode a 
In the case of a sheet beam, this qua si-electrode (the sheet beam) may be 
pla cod closer to the source without voltage breakdown since extracted ions 
are free to pass completely through this extraction beam. Another advan- 
tage of using a sheet of electrons as the extraction electrode is that there 
is no longer any erosion which might occur due to the sputtering of some 
of the ions onto a solid electrode. (Again, the ions pass freely through 
the sheet electron beam.) Charge neutralization represents another way of 
increasing the current density. If particles of sign opposite to those 
being extracted (for instance, electrons) are injected in the extraction 
region, neutralization of the resulting space charge is affected which, 
in terms of the Langmuir-Child law , effective^ increases k, the quality 
factor. This unwanted charge may be separated from the desired ionic 
charge later by means of magnetic fields (taking advantage of the great 
disparity between electron and ion cyclotron radii) as is done in isotope 
separators. Still another way of increasing the current density, J, is 
suggested by Danilov (11, 12). In his papers, Danilov states that by 
properly shaping the magnetic fields in the extracting region, it is 
theoretically possible to again increase k. This is a slightly different 
kind of space charge neutralization than that mentioned above in that in 
this case, neutralization results from magnetically induced Hall currents 
which will flow in the extraction region. These are three ways which will 
be investigated in the present research in order to try to maximize the 
output of the ion source. 

The source of ions which will be utilized in the proposed research is the 
magnetron source. As was mentioned earlier, Masek (7) found that the 
inverted magnetron yielded a more uniform ion density than the normal 
con f iguraticn . A great advantage in using this kind of a source as 
opposed to a contact cesium-tungsten device is that corrosion is greatly 
reduced in using the magnetron. In his contact source, Meyer (2,3,4) 
reported that corrosion of the ionizer was indeed a major problem. This 
problem will be greatly reduced by use of the magnetron as the possible 
source of the plasma from which ions will be extracted. 

A First Theoretical Model 

As is obvious, the complete field solution of sheet beam of electrons 
passing over a plasma, formed in an inverted magnetron, is extremely diffi- 
cult to obtain. Therefore, the first attempt at understanding the physics 
involved in this problem is contained in a hydrodynamic model for the 
plasma. The plasma is, therefore, approximated to be a highly conducting 
inviscid incompressible irrotational fluid with no free charge, which 
satisfies the Na vier-Stokes and conservation! of mass equations. Also the 
sheet beam is approximated by an electrode of potential, -Vo> with respect 
to the plasma. This seems to be a reasonable assumption since the poten- 
tial of the electron sheath would be approximately that of the tungsten 
cathode which is part of the electron gun. In addition, since the operat- 
ing pressure is approximately 10*"^ torr, the mean free path of these elec- 
trons is large and, therefore, they do not interact much with each other. 
The model of this system is shown in Figure 1. 
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Figure 1» 



The dispersion relation for 
tion technique as described 
perturbations of the form 



this system is found by the use of a perturba- 
by Chandrasekhar (13) and Melcher (14) * Small 




are assumed about an equilibrium which occurs when the fluid surface is 
at X = 0 and the electric field is ^e> . The solution for the 

perturbation fields which exist when the fluid interface deforms slightly, 
and for the fluid velocity and pressure is derived explicitly in Melcher *s 
text (14). The resulting dispersion relation is 



o>^ = 






where 






(3) 



te - \ k'/ + 



2 ^ 

X 



and demonstrates exchange of stabilities 
This means simply that for a given real 



the wave number in the x-direction, 
as defined by Chandrasekhar (13) . 
value cf k, instability occurs for 



O and the instability is absolute; that is, it is a growing 
exponential in time, not a growing oscillation. It is important to note 
that by definition if even just one mode (i.e., one value of k) leads to 
an instability for a given disturbance, the system is said to be unstable. 
Onset of instability occurs for O in which case equation (3) 

becomes 
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( 4 ) 



where x = kdo It is interesting to investigate this equation in its two 
basic limitSo They are the short and long wave limits in which case the 
value of voltage, Vq) necessary for instability, is given as follows: 



short wave limit 



I 



long wave limit 







where ^ is the mass density, T is the surface tension and is the per- 
mittivity of free space., The resulting dispersion relations for each 
limit are shown in Figure 2 as a function of Vq • 




In an attempt to corroborate this theory the following experiment was 
performed. A liquid metal, the eutectic f allium-indium alloy*, which is a 
liquid at room temperature, was stressed by a perpendicular field much as 
in Figure 1. As the voltage was increased some of the metal was pulled 
toward the upper plate. This initial disturbance quickly grew into a 
sharply pointed spout. Instability was assumed to occur when this spout 
occurred (just before arcing to the upper electrode) since the initial 
disturbances, when the metal was just pulled toward the upper electrode, 
would eventually die out with a long time constant. Voltages for the onset 
of instability are plotted in Figure 3 with the theoretical result from 
equation (5b). The instability is found to be of the long wave limit since 
the surface tension of this liquid is so great. Melcher (15) also finds 

* 

This metal had the following parameters in mks units: 
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^ = 6700, T = .63. 







this type of result because the factor G = Pgd2 , which is the ratio of the 

T 

energy density stored in the gravitational field to the energy stored by 
surface tension, is less than 3, (It is 2.6 in the case at hand<») In 
liquids for which this is true, dispersion relations of the kind found in 
Figure (2b) result . 



There are several reasons for the disparity between the theoretical 
and experimental results as shown in Figure 3, It was assumed initially 
that equilibrium occurred when the fluid interface was flat and stationary. 
It was very difficult, however, to maintain this kind of a boundary as the 
point of instability was approached so that one of the initial assumptions 
was violated. Also, the surface of the gallium-indium alloy was difficult 
to keep "clean. Impurities of some kind, perhaps oxides, formed on the 
surface and were large enough at times to distort the surface of fluid. 
Thus, these distortions resulted in areas of increased field strength. 

All attempts at purification failed. As predicted by theory, the first 
wavelengths to become unstable were long waves. The experiment seemed to 
predict this since instabilities occurred at only one point on the inter- 
face and not at regular intervals which would suggest a shorter wavelength. 
Since the experimental apparatus was only 2.5 inches square (due to a lack 
of sufficient gallium-indium alloy for a larger device) , the effect of the 
boundaries on the fluid interface would probably be great. In addition, 
the experiment was perfonned in vacuum and the apparatus viabrated slightly 
along with the vacuum equipment. These viabrations never did allow the 
surface to be flat and the separation between the initially slightly 
rippled Interface and the upper electrode was effectively less than 
measured. All of these sources of experimental error would tend to cause 
the interface to become unstable before the theoretical value. 



This paper was supported in part by a grant from the U.S. Air Force, 
AFOSR-68-1508. 
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THERMONUCLEAR REACTIONS INITIATED IN PARTICLE- PARTICLE COLLISIONS 



C. D. Hendricks, A. W. Dipert and D. Weidenfeld 
Charp^ed Particle Research Laboratory 
Department of Electrical Engineering 
University of Illinois 



I nt roduct ion 

The possible production of a controlled thermonuclear reaction 
initiated from a solid-solid interaction at extremely high velocity and 
laser beam irradiation is currently being investigated. The production of 
a thermonuclear reaction requires that the interacting particles have a 
temperature equivalent of about 10 KEV per nucleon. Calculations for a 
realistic accelerating potential will show that there is an upper limit 
to particle size for a fusion reaction initiated only by high velocity 
particles. If large size particles are used, or a lower accelerating 
potential is applied, the additional energy required for a fusion reaction 
may be supplied by irradiating the particles at the moment of impact with 
a high power laser beam. Preliminary experiments will be done with a 150 
KV power supply. Since the supply potential is not nearly sufficient for 
obtaining fusion velocities for any reasonable size particles, a high 
power laser will be used to supply the addit iona 1 required fusion energy. 

Basic Theory 

If a charged particle is to ignite a thermonuclear reaction, the 
atoms involved in the collision must have an energy of at least 10 KEV. If 
a particle (Lithium Deuteride is to be considered) is acceleraged through 
a voltage of 10,000 volts, it will have sufficient energy for ignition of 
a nuclear reaction provided that the particle is 100 per cent ionized. 

Since this is impossible for a solid particle composed of more than a few 
atoms, we will calculate a lower limit of ionization given a practical 
accelerating voltage. Many electrostatic generators are capable of pro- 
ducing a potential of at least 10 million volts, and if we set this value 
as the upper limit for accelerating voltage, then at least 1 out of a 1,000 
atoms must be ionized to achieve sufficient fusion energy. 

The required cha rge-to-mass ratio for a particle with a fraction 
(f) of its constituent atoms ionized is given by the expression 



Q ^ fZN^qp 
M A 

Z = number of protons in a molecule 
Nq = 6.02 X 10^^ molecules/kg-mole 
A = molecular weight 
q^ = 1.6 X 10“"^^ coulombs 

Then for an f value of 0.001 for a LiD particle, the charge-to-mass ratio 
is approximately 2 x 10"^ coulombs/kilogram. 

Next, the electrostatic field E^ at the surface of a spherical 
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pnrticle is computed for Lithium Deuteride for a set of radii* for a 

surface of a sphere is given by 

- Q 
4 n r^ Eq 

Q = total charge = Nqo 
N = Number of ions/particle 

Also using the equation for number of atoms per kilogram-mole and the 
equation for the mass of the particle 



N = ^ 
A 



m 



4 3 

m =-= n 



Combining the above equations, an expression for Eg in terms of the 
particle size (radius) and physical properties of the particle can be 
obtained. 



p = Npqn 

s 3 

o 




For fractional ionization, Eg’ = Egf » and using the above equa- 
tion and f = 0.001, the surface field for LiD is computed for a range of 
particle sizes. 
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Z = 2 

P = .883 g/cc LiD 
A Cy 9 



Using a value 10^ volts/meter given by most references as an 
upper limit to the electric field at the surface due to field emission 
effects, we require a particle size of about 0.01 microns. 

Development of a Reaction System 

The success of a system to produce a reaction depends on the 
development of a method for obtaining a uniform rate of charged solid 
particles with a consistent charge-to-mass ratio. Current research efforts 
are investigating particle charging methods, charge-to-mass ratio selection 
schemes, and particle containment and ejection systems. 

The quadrupole electric field is investigated as a practical 
method of particle selection and containment. In the device, the particles 
are held in dynamic equilibrium by alternating electric fields. The 
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electrode configuration of the device gives sinusodally time varying 
forces whose strengths are proportional to the distance from a central 
origin. The governing differential equation of particle motion is a special 
case of the Mathieu differential equation. The solution of the Mathieu 
equation foi* stable motion is dependent on the charge-to-mass ratio of the 
particle and the applied driving frequency and voltage of the device. By 
selecting an operating voltage and frequency, a range of stable charge-to- 
niass val ues can be calculated. 

Experimental Apparatus 

A small experimental electric quadrupole chamber was machined from 
brass (see diagram). For initial experiments, 60 cycle line frequency was 
used. A high voltage transformer connected to an auto-transformer gave a 
variable voltage between 0 and 5,000 volts. The chamber was mounted over 
an electric particle injector. The particle injector is essentially two 
parallel plates spaced about 1 cm apart. The upper plate has a small 
aperture for particles to escape. A high voltage is applied between the 
plates and the electric field charges and attracts the particles to the 
upper plate. A fraction of them pass through the aperture and enter 
the quadrupole chamber. 

Particles with the correct initial velocity and charge-to-mass 
ratio should then be retained within the chamber. In the actual experiment, 
after a few seconds of operation, between 1 and 5 particles were observed 
to be suspended. By varying the driving potential Cor a specific value 
of charging voltage, an optimum value could be found that would give a high 
suspension rate. By switching a high voltage on the upper plate of the 
suspension chamber, the suspended particles could be ejected out the small 
aperture in the upper plate of the quadrupole. 

The experiment showed that a relatively simple system could be 
constructed to select and eject charged particles. Refinements of the 
system could be made to increase the ejection rate, narrow the charge-to- 
mass selection, and increase the charge per particle. 

Conclusion 

The current exploratory experiments have completed the initial 
design of a solid particle injection system capable of producing a uniform 
rate of charged particles with a consistent charge-to-mass ratio. The 
next area of investigation is to refine the system to obtain a higher 
charge-to-mass ratio and a higher particle ejection rate. 

Future charging investigations would use electron beams to increase 
the charge on a suspended particle. Circuitry would be developed to vary 
the driving frequency of the device as the charge increased. Studies of 
high power laser (in the gigi-watt range) characteristics will be conducted. 

A suitable laser system will then be purchased or constructed. A 150 KV 
power supply will be used for acceleration of particles. Timing techniques 
have been developed to synchronize the laser pulse with particle collision. 



This paper was supported in part by a grant from the U.S. Air Force, 
AFOSR-68-1508. 



25 



DIAGRAM OF THE PARTICLE 
SELECTION-EJECTION DEVICE 
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INTERACTION OF LASER PRODUCED PLASMAS 



WITH A MAGNETIC FIELD 

M. J. Lubin and H. S. Dunn 
University of Rochester, Rochester, N. Y. 

For the past three years the plasma physics group at the 
University of Rochester has been investigating the interaction 
of laser produced plasmas with strong magnetic fields in both the 
collisionless and collisional regimes. Use of the energy of 
giant pulsed lasers to create a plasma of known and measurable 
properties, devoid of externally applied driving currents, has 
proven to be a unique experimental tool for the study of plasma- 
magnetic field interactions and nonlinear wave propagation. Here 
we report on two phases of our work, the electromagnetic radiation 
from a laser produced plasma expanding into a uniform magnetic 
field and the two-stream instability (and possible collisionless 
shock wave) created by the expansion of a laser produced plasma 
into a quiescent background plasma. In both projects, experimental 
observations have been supported by theoretical studies. 

1. Electromagnetic Radiation from £ Laser Produced Plasma 
Expa nding into ^ Uniform Magnetic Field (1). Experimental 
observations (2) of the radiation field resulting from the 
expansion of a laser plasma produced from the vaporization and 
heating of a droplet of lithium hydride in the presence of a 
uniform magnetic field suggest the existence of strong azimuthal 
surface currents. A typical oscilloscope trace of the radiated 
power is shown in figures 1 and 5. The radiation measurements. 
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all made in a direction normal to the magnetic field as shown in 
figure 2, were carried out using a submillimeter detector with 
measurable wavelength response from 70 microns to 8 mm. combined 
with a rise time of less than 4 nanoseconds. 

The result of a computer simulation of a representative 
experimental situation where the magnetic field is sufficiently 
weak that the expanding plasma remains spherical during the time 
interval of interest is shown in figure 3. To calculate the 
electromagnetic radiation signature of the expanding plasma, it is 
necessary to consider three separate effects, namely the dipole- 
like radiation resulting from the induced azimuthal currents, 
the surface cyclotron radiation and volume cyclotron radiation. 

The dipole radiation field is obtained by solving Maxwell's 
equations both inside and outside the expanding plasma sphere and 
matching the two solutions at the plasma-vacuum interface, r=R(t). 
The external solution is easily obtained and the far-field power 
radiated into a solid angle c\il is 






"f). IF") ^ ^ 



( 1 ) 



where F is an arbitrary function of the retarded argument t-r/c 
and the primes on F denote differentiation with respect to the 
argument. For the internal solution, we use directly the results 
given by Bernstein and Fader (3) to deduce, after many approximations, 
that 
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( 2 ) 



- ' 




for values of R '^/2, the length being given by 




(3) 



where the subscript "a" refers to the asymptotic phase of the expan- 
sion. Inserting this value of F into the far field approximation 
eq. (1) gives the radiated power. 

As the plasma temperature drops during the expansion, penetra- 
tion of the magnetic field into the plasma volume begins to occur. 
During this time, a certain number of electrons escaping from 
the plasma will be reflected back into the plasma volume through 
helical trajectories. Multiplication of the average intensity of 
radiation for a single electron in cyclotron motion (4) by the 
number of electrons escaping from the plasma volume gives, after a 
detailed calculation, the total intensity radiated into a solid 
angle l\A. , 



where it has been assumed that observation is normal to the magnetic 



After the magnetic field has diffused into the volume, we 
would expect volume cyclotron radiation to occur. Integration 




(4) 



field and 
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of volume cyclotron radiation from an elementary volume dT of 
plasma (5) over the plasma sphere gives 



e\l 

\3 ^ 




( 5 ) 



for observation normal to the magnetic field. 

The total contributions of the dipole, surface and volume 
cyclotron radiation are plotted in Figure 4. The conductivity 
in these calculations has been the scalar conductivity associated 
with the free expansion. During the peak of the radiation signal, 
the r correction is significant since \ t is of 

the order of 10 or higher, thus leading to a smaller value of "S~ 
which should appear in equation (3). This effect is also offered 
as a tentative explanation for the trend of the measured radiation 
intensity in figure 5 toward a linear dependence for higher values 
of magnetic field. The radiation intensity illustrated in 
figure 4 was obtained by adjusting the relative contributions of 
the dipole and surface cyclotron radiation to fit the general 
shape of the experimentally measured form. A detailed fit is 
difficult. This is understandable in view of the crudeness of 
the model used to calculate the radiation history, the limited 
experimental data available, and the lack of knowledge of the 
surface number density required to calculate the surface cyclotron 
radiation. A more meaningful problem would involve the study of 
the interaction with uniform magnetic field when spherical symmetry 



30 



of the expanding plasma is not assumed. This is presently being 
pursued. 

2. Two - Stream Instability in Laser Produced Plasmas (6). 

Freely expanding laser produced have large directed energy 
to thermal energy ratios. The expansion of such a plasma into a 
quiescent background plasma is being studied both experimentally 
and theoretically to investigate the properties of two-stream 
instabilities and possible existence of collisionless shock 
waves (7,8). 

Our experiment utilizes a cesium, Q-machine plasma, as a 
background for the expansion of an aluminum plasma formed by 
focusing a 40 MW, 20 nsec, wide ruby laser pulse onto a 20 micron 
thick aluminum wire. The hot plasma produced from the vaporization 
of the aluminum wire rapidly cools as the thermal energy absorbed 
from the incident radiation goes into directed energy of expansion. 
After a few nanoseconds the aluminum plasma is expanding with its 
asymptotic expansion velocity. The Q-machine used in this experi- 
ment is a single ended, pulsed emitter type, similar to one built 
at Columbia University (9). The laser beam is focused onto the 
wire which is located at the center of the cesium plasma column. 
Alignment is done by means of an auxilliary gas laser. The 
expansion of the laser plasma is monitored by floating double 
probes located along the axis of the Q-machine. 

The experiments are conducted in a low axial field (B=100 gauss) 
in order to minimize the effects of the magnetic field. The results 
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of the experiment, summarized in figure 6, show four typical probe 
traces taken at different background conditions. A definite.., 
decrease of the peak current was noted as the background density 
was increased. Apparently the presence of the cesium plasma caused 
a number of the free streaming aluminum ions to be either deflected 
or slowed down. 

The mean free path for coulomb scattering is given by: 



particles and n is the density of scatterers. Since the measured 



particle scattering cannot be explained on the basis of simple 
Coulomb scattering. 

An analysis was undertaken to find out whether the experimental 
results could be explained on the basis of a two-stream instability. 
Following Stringer (10) the dispersion relation for a many component, 
one-dimensional, homogeneous plasma may be written as: 




\\ \v^ A- 

where W is the energy of relative motion of the two interacting 



aluminum expansion velocity is about 5 x lO^cm/sec., the cross 
section for aluminum-cesium coulomb collisions is very small 
yielding mean free paths on the order of 10 cm. Thus the observed 




where AL*' 
.1 




4lT\V, /OV (plasma frequency of each component) 

J J J 

sound speed of each component 



V.’, 



drift velocity of each component 



Xl2r) ■= i (jjU) 
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where is the error function of complex argument. 

In applying this relation to our experiment the following assumptions 



were made : 

1. The electron and ion-temperatures of the cesium plasma 
are equal. 

2. The electron and ions in the laser plasma are expanding 
at the same velocity. 

3. Since the temperature of the cesium plasma is T=.25ev., 
the electron thermal velocity of the cesium plasma was 
taken to be five times greater than the expansion velocity 
of the aluminum plasma. 

The dispersion relation was solved for various electron and ion 

temperatures in the laser plasma and for several density ratios of 

cesium plasma to laser plasma. For a given density ratio, plots 

of the instability regions and various growth rates were made as 

functions of two parameters, the electron, ion temperature ratio- 

Tg/T^, and the ratio of directed velocity to electron thermal 

speed for the expanding plasma, V/Ce. A typical instability 

plot is shown in figure 7 for a density ratio (aluminum density/ 

cesium density) Nj^/N 2 = 10. For a streaming velocity of 4 x 10 

cm/sec. , a representative growth rate is ^ • — ’ 5 x 10“^_[1_ pe 

which gives interaction distances of approximately 10 cm. for a 

11 -3 

streaming density of 4 x 10 cm . Measured interaction distances 

are from 6 to 10 cm. this giving reasonable agreement. 

From the measured probe currents, we estimate the density of 

11 3 

the expanding plasma to be no more than about 5 x 10 /cm thus 
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giving typical density ratios of 10 to 100. The growth rate of 
the instability versus the density ratio is plotted in figure 8 
for the particular case = 10^ and given V/Cg. As a matter 

of interest we plotted the number of particles collected in the 
main peak of the probe traces versus background density, this plot 
is shown in figure 9. In this plot N refers to the particles 
collected with background cesium and Nq to the number of particles 
collected with no background plasma, In figure 9 we see that the 
number of particles lost roughly follows an exp dependence 

characteristic of a two-stream instability rather than the exp(^-^bn) 
dependence to be expected for coulomb scattering, again re- 
affirming the fact that coulomb scattering of the ions is not 
responsible for the observed particle loss. 

The experiment then leads us to the following conclusions: 

1. Laser produced plasmas can be two-stream unstable even 
when expanding into much lower density plasmas. 

2. Collisionless shock waves will be difficult to produce 
using laser plasmas unless much higher density background 
plasmas are used than were used in this experiment. In 
the case of higher density background plasmas snow plowing 
can be expected to play an important role. 

The effect of a strong axial B-field in the interaction of 
the two plasmas is presently being studied. 
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FIGURE CAPTIONS 



1. Oscilloscope Trace of Far Infrared Radiation Emitted from 
Laser Plasma Expanding in a Uniform Magnetic Field 

2. Geometry of Expanding Plasma for Dipole Calculations 

3. Computed Time-History of Expanding Laser Plasma 

4. Computed Radiation History 

5. Measured Values of Far Infrared Radiation Intensity in a 
Uniform Magnetic Field 

6. Typical probe signals due to laser produced plasma 
expanding against background Cesium Plasmas. Sensit- 
ivity is 10 ma/division, .5 sec/division, 

(a) no background plasma, (b) background density of 
4'x 10^/cm^, (c) background density of 1.2 x 10^®/cm^, 

(d) background density of 2.7 x 10^^/cm'^. 

7. Stability conditions for aluminum streaming into a back- 
ground cesium plasma, density ratio (aluminum/cesium) = 10. 
The parameters V/Ce and Ti/Tg and the plasma frequency, 
refer to values for the streaming aluminum plasma, 

is the growth rate of the fastest growing waves. 

8. The maximum growth rate for Tp/T^ = 2.5 x 10 plotted 
as a function of the ratio N 2 /N^ = (cesium/density/ 
aluminum density). 

9. Measured particle loss versus background cesium plasma 
density. 
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LASER CREATED PLASMA AND COLLI5IONLESS SHOCK WAVE PHENOMENA 

by 

F. Schwirzke and A. W. Cooper 
Naval Postgraduate School, Monterey, California 

1 . Introduction 

Collisionless shock waves are an important mechanism in heating plasmas 
to temperatures and densities of thermonuclear interest. Laser production of 
plasma is now being applied to the problem of investigating collisionless 
shocks in a laboratory plasma. A steady state collisionless plasma is formed 
in a hollow cathode arc facility. Collisionless shocks are produced by the 
rapid expansion of a hot plasma obtained by focusing a giant laser pulse on the 
tip of a small target. The expanding laser plasma acts as a very fast piston to 
produce the collisionless shock. The adjustable experimental parameters 
include background plasma density, external magnetic field, ion mass, Alfven 
velocity, electron temperature, and the velocity of the expanding plasma. By 
thus varying the magnetic Mach number the structure and damping mechanisms 
of the collisionless shock can be studied in an effort to deduce the collisionless 
mechanisms of energy dissipation and its effectiveness in thermalizing the 
ordered wave motion. 

2 . Dynamics of a Laser Created Plasma Expanding in a Magnetic Field 

The two problems of creating and confining a high temperature plasma in a 
magnetic field configuration are usually closely linked together and a laser 
produced plasma is no exception. Theoretical calculations and experimental 
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results^^^®^ show that energy from a giant laser pulse can be absorbed by a small 
target within a few nonoseconds, thus creating a high temperature plasma. 
Obviously the next step is then to investigate the interaction of the expanding 
plasma with a magnetic field. This interaction is of special importance in the 
controlled thermonuclear research program where the purpose is to confine a 
thermalized hot plasma in a stable magnetic field configuration. In the experi- 
ments to be discussed^^®^ the plasma was produced by impact on a 30 micron 
glass fiber of a 20 nsec, 150 MW pulse from a neodynium glass laser with no 
background plasma. About one joule of energy was absorbed in a plasma of 10^® 
atoms, giving an expansion velocity of about 10^ to 10® cm/sec, and initially 
a plasma-to-magnetic pressure ratio of order 10^. 

a. The expansion of the B 10^ plasma can be described by hydrodynamic 

theory. The temperature decreases rapidly according to the form T = T (r /r)^, 

r o o 

where T and r are the temperature and radius of the plasma when it ceases to 
o o 

absorb laser energy. Most of the energy content of the plasma appears as 
kinetic energy of expansion. The rapid decrease of the plasma temperature implies 
that an asymptotic plasma expansion velocity is reached in a rather short time 
while B >> 1. This explains the observed asymptotic expansion velocity being 
the same with and without a magnetic field. 

b. The expansion of a plasma in the presence of a magnetic field has been 

treated theoretically, and an oscillation or bounce of the plasma against the field 
predicted at the condition of B Experimental observations are in approx- 

imated agreement with this result^^°^. An expanding plasma of infinite conductivity 



47 



would repeatedly bounce off the compressed B where p ^ 1. With finite resistivity, 
diffusion across B will also occur. Hence, the bounce radius and the character- 
istic width for the pressure gradient in the wave front are determined by the 

effective resistivity. This means a laser produced plasma which expands across i 

1 

a magnetic field also reveals extensive new possibilities for the investigation of 

instabilities and turbulence in the wave front^^^’^^\ 

c. Resistive effects associated with surface currents in the wave front form 

the most effective means of rethermalization of the directed kinetic energy of 

expansion. In addition, if the growth rate becomes greater than the electron 

Coulomb collision frequency u / instabilities may develop in the front and also 

e 

contribute very effectively to the energy dissipation process. 

The two-stream instability will start where the current-associated mean 

ordered velocity of the electrons in the wave front becomes larger than the mean 

thermal velocity. The current is in the direction perpendicular to B and to the 

normal of the wave front. The two-stream instability has a large growth rate 

y cc' (m/M,)^^^ and may develop if p < y , where to is the plasma frequency 
pe * e pe 

and m and are the electron and ion mass. The growth time of the two-stream 
instability is rather small in comparison with the propagation time of the wave 
which is in the order of 



T ^ 10 

Here V. is the Alfvdn velocity, d 10 c/to is of the order of the characteristic 
width of the front when the resistivity of the plasma is determined by the presence 
of the instability, c is the velocity of light. If yT>> 1 a steady state will exist 
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within the wave front and the onset of the instability can be studied locally by 
considering the perturbed motion of the ions and electrons by means of the equations 
of the two-fluid theory. 

Considering a stationary wavefront which is common to shock theory, 

(see figure), propagating in the r-direction across the magnetic field, the electric 
current is perpendicular to B and r in the 0 direction and is given by Maxwell's 
equation 

j = -enV^Q = - (c/4tt )dB/dt. 

Substituting (x^ = Anne^/m and K = | (l/B)dB/dr) i where K is inversely 

pe 

proportional to the width of the front d, the criterion for the onset of the two- 
steam instability, ^ kT^/m , 

becomes: d < (ttc/co^^) (2/0)^^^ (1) 

where d = (8tt nkT^/B^) . Criterion (1) will be satisfied at a certain radius since 
n and decrease rapidly during the expansion of the laser produced plasma. In 
agreement with these discussions the measurements indicate that the two- 
stream instability enhances the resistivity in the wave front by more than two 
orders of magnitude and thus determines the width of the wave front of about 
d ISc/to^^. The discrepancy between the values of T obtained from the magnetic 
probe signal and from the effective conductivity as estimated from the width of the 
front also indicates that a collisionless mechanism is contributing to the energy 
dissipation. 
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RUCTURE of the wave front (SCHEMATIC) 



3 . Conclusions 



Theory and experiment have shown that a laser produced plasma can be used 
to study turbulent energy dissipation mechanisms and collisionless shock wave 
phenomena. However, experimentation requires higher power levels, and 
therefore the laser system is now being expanded to about 1 GW by addition of 
an amplifier stage. 
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STREAMING PLASMA INTERACTION WITH 
LONGITUDINAL MAGNETIC FIELDS 



I 



J. G. Siambis 

Carnegie-Mellon University 
Pittsburgh, Pennsylvania 15213 

It is well known that plasmas produced in theta pinches satisfy 
thermonuclear power requirements in two respects, density and ion tem- 
perature but not confinement time. In this work the transfer and flow 
of a plasma from a 5 KJ conical theta pinch gun is studied as it enters 
a rising longitudinal magnetic field. The maximum value of the single 
mirror magnetic field can be raised to 7 KG. The following diagnostics 
have been utilized, 70 GH microwave interferometer, STL image converter 
camera, spectroscopy, diamagnetic loops, photomultipliers and probes. 
Measurements have been carried out at the location of maximum magnetic 
field gradient. We find that the electron density rises from 10^^ per 
cc for zero field to 7 x 10^^ per cc for 3 KG at the measurement loca- 
tion. The electron temperature rises at first from 6.7 ev to a maximum 
of 12 ev at 1 KG and then falls to a value of 2 ev for the maximum 
field of 3 KG, while the ion temperature as obtained from the diamagnetic 
loop measurement increases throughout this range. The front velocity 
of the plasmoid is found by probe measurements to be 5 x 10 cm/sec. 

One of the objectives of these measurements is to establish the occur- 
ance of a shock wave in the transition from supersonic to subsonic 
flow in the plasmoid. (see Figure 1 ) Further modifications in 
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the system will incorporate a guide field between the gun and the mirror 
field and thus minimize 'the interaction with the tube walls. 

The theoretical work is concerned with the propagation and 
stability of drift waves in realistic magnetic field configurations. 
Results have been published^ which elucidate the effect of shear on the 
stability of drift waves in stellarator and levitron type geometries. 
Further work in progress takes into account the variation of magnetic 
field along the field lines and drift-Alfven waves are also considered. 




Figure 1. Image converter camera pictures of plasma flow from plasma 

gun (right) to converging longitudinal magnetic field (left). 
The time is shown in microseconds with each picture. The 
converging magnetic field to the left has a value of 3 KG, 
and is the location of maximum magnetic field gradient. 



^J. G. Siambis, Phys. Fluids 12, 1971 (1969) 
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Vorticity and Neutrons in the Plasma Focus 
W. H. Bostick, L. Grunberger, V. Nardi, and W. Prior 
Stevens Institute of Technology, Hoboken, N.J., U. S. A. 

Image converter photographs (5-nanosecond duration) of the current sheath 
in a plasma coaxial accelerator operated with at 8 torr show that the current 
sheath is corrugated and that there are imbedded in this sheath dense plasma 
filaments (see figures 1 and 2). Five-nanosecond time-and space-resolved spec- 
troscopy shows that the light from these filamentary structures appears 

primarily as a continuum. The degree of broadening of the H3 in the filaments 

18 3 

places the 'electron density there > 10 /cm . The stream velocity ^ (relative 
to the current sheet) for mass flow ahead of and on the current sheath is of 

the same order of magnitude of the current sheet velocity in the laboratory 

7 1 

frame (yO.SxlO cm/sec). Evidence has been adduced to suggest that these 

filaments, which are produced in pairs, coincide in space with regions where 

the magnitude |o)| of the vorticity co = 1/2 curl jv assumes very high values. 

We refer to these filaments as plasma vortex filaments. Mass flow along the 

filaments from center conductor to outer conductor has been observed. Magnetic 

fields parallel and anti-parallel to this flow along the filaments have been 

measured. These magnetic fields occupy regions of space comparable in extent 

to the size of the filaments and are associated with the filaments. Magnetic 

probe measurements indicate that a large portion of the total current in the 

plasma coaxial accelerator is carried in the current sheath, and that the current 

sheath is about 2.5 mm. thick as compared with the filament diameter, 
which is about 1mm. The corrugated shape of the current sheath (figures 1 and 2) 
is an indication that a large fraction of the current carried by the sheath is 
actually carried by the filaments. 

Schlieren photographs taken with a giant pulsed laser with an exposure 
time of about 25 nano-seconds show essentially the same features as the image 
converter camera shows and place the electron density in the filaments at a 
value of about 10 times the background density, or a filament density of 
3xl0^®ciri^. 

When the spacing between the vortex filaments is small enough, they 
annihilate in pairs, and at the plasma focus (figure 3) near the axis a 
general region of annihilation (enhanced illumination) starts at the end 
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of the center conductor and moves away from the center conductor along the axis 

7 

at a speed of 1.5 x 10 cm/sec. 

Two main neutron pulses occur: the first during the vortex annihilation 
along the axis while the bright region moves from the center conductor to 
the crown of the ’’umbrella”. The second main neutron pulse is associated 
in time with the vortex annihilation occurring in the crown of the umbrella 
as the vortices ’’burn” pair by pair as a shock wave moves radially outward 
in the crown. This observed vortex annihilation is the authors believe, the 
laboratory scale version of the solar flare phenomenon. With this relatively 

g 

small scale apparatus the neutron yield per pulse is about 5 x 10 neutrons . 

6 7 

This yield is much less than that obtained by other experimenters ’ who use 
much higher energy storage. 

On the basis of the experimental observations it seems reasonable to con- 
clude that the observed ’’optical front” represents very well the geometric 
configuration "the region (of thickness 1mm) of sharp transition for mass 
density p pressure (isotropic) P, fluid velocity and magnetic field 

Specifically, by considering as usual the shock front in the limit d ^ 0, 
i.e. as a surface of discontinuity Z (say) for these quantities, the most 
appropriate location of Z is obtained by identifying Z with the foremost face 
of the optical observed layer. The filaments which appear embedded in the 
optical front, i.e. immediately behind Z, can be related closely theoretically, 
as we show here, to the shape of Z. 

Since the equations for the magnetic field and for the vorticity become the 

2 

same under some circumstances and even a relationship of proportionality 

between magnetic field and vorticity can be deduced for finite-conductivity com- 
3 

pressible plasma it is suggested that a similar relation can exist also for 
different or more general conditions and in particular for our flow-carrying 
filaments and the magnetic fields parallel or anti-parallel to the filaments. 

The vorticity, V x v^, behind the shock front (suffix o denotes conditions 

ahead of the shock, suffix 1 conditions behind it) can be evaluated in terms of 

the two ’’principal curvatures” of Z ; in terms of two shock strength parameters 
P ]_ 

C = , a(H); and in terms of v^ , . 
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A convenient form of the Newton’s equation for stationary conditions (3^=0) 
and infinite conductivity ^ ii - 'the electric field in the plasma) is 

(1) U X li) X = - -ij I 5 I vj - X (V X H^) 

The system of orthogonal curvilinear coordinates with unit vectors , x, 

a D ^ 

tangential to the ’’curvature lines” of E and n_ normal to E, is introduced 
(A A . . 

’ “XI then indicate respectively the tangential and normal components of 

a vector _A with respect to E) . 

Since [(V x v ) x v 1 = (V x v ) x v^ and 

— —1 —lx — ~1 —Ijn 



(2) V X V = - [(V X V )xv ] ^ T7 
~ —1 — —1 —l,n ^ 



l,n 



, the tangential component Xx_^ + of 



the right hand side of (1) is first estimated in order to obtain V x v^. 

All quantities p_ ,v , H . P can be expressed in terms of p ,v , H 

1 — ^1 — 1 1 0 — o — o o^ 

expressed by the jump (shock) conditions on E (deHoffman, Teller relations ). 

The leading contributions to X,y (the differential operators are transformed to 

curvilinear coordinates and ’’strong shock” conditions are assumed, P "^^o , 

o 

enthalpy i^ ^ 0) have the form 

where ^ and ^ are the ’’principal curvatures” of E (at a point x the curvatures 
K K, — 

are taken positive when they are convex to the incoming flow) and 

Hi a ^ 

ot r - — 1— = - — if ^ > 1. By substituting in (2) we have 
“o,a “o,b 

1 - T — — + X — 

(3) u) = — V X V = — a 2v R, ^ 2v R 

— 2 — —1 o,n D o ,n a 



If X is the vector with projection on the orthogonal-to-z plane, which is(approx- 



— a 



imately) collinear to the electrode radius R then, we have by experimental 



observation (see figures 1 and 2), that 



. '"b, 



up to and during the first 



stage of the ’’rolling off” of the optical front onto the face of the center con- 
ductor. Here the filaments which are perpendicular to the azimuthal magnetic 
field predominate. Later at the final stage of collapse increases and 

decreases and it is now possible to see filaments in both t^e azimuthal and radial 
axial directions (figure 3, note, ’’spider web” type of filamentary pattern). 
Correspondingly, during the earlier period, o) , which is collinear with the radial 

Si 
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♦ 1' P' ^ J i I . 



t:*e 




axial fiiamont®, dominant, but r , •M*.ri 

azimutyial filair.ont , an-' to_ ar . . 



The algebraic 3igtiS oi trie w componerut r *‘gnj t R_ , k 

ri C 

(and also on the functions F,C). ?ot', : ign* ai*- ;lt^ir^^y i.ehind with 

a relative frequency ipparenliy In agricn^nt v/iiii cne eirfti.riL.al o.i^ervations . 

In r uirmary , e-] la’^lcn (w , ’^r* ^ r; n c^rr^ >1: ^n* ±iiV^ ana '-j^-- In ^as ../na-nlc 

r 

shocks )jhow .• that th(. exL^t^n,^ tf vortlcity pjr-l.,ei t t: o ciFervz-. liga- 
ments can be pi’edic ted as a Trc'^t -it’riren of tii* c r ’^■.itore of tnc 

ob s e r ve d shock i r ‘o n t . 



In add’tion t-^ he MHI3 sol«itfrn ^'ne of t:i'* aut ir- ^ .♦!€ ♦ .t .inc*.- a /cuu+Aon 
to the Vlasov equation which givet tt'i* < alr. j* f i^cjr«fi'. .. wr?f. lHc cxpe^ .ed 
r'^"'s flow tr iL'l irc ma lt(»A^nctic •!* . iiv . 
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Figure Captions 



Figs. I5 2, and 3 are 5-nanosecond exposure IC photographs of the PCA operated 
with 8 mm of deuterium. The time at which each picture was taken is indicated in 
terms of t = 0 coinciding with the first peak in the ^ signal. Figures 1 and 3 ; 

are taken along the electrode axis z. Figure 2 is a profile. Note that the "scallops 
appearance of the current sheath in figures 1 and 2 is evidence that even at this 
early time the vortex filaments are carrying a large fraction of the total current 
in the sheath. The PCA has inner and outer conductor diameters of 3.4 and 10 cm, 
and an inner conductor length of 13.6 cm. It is driven by a 45 yf capacitor bank 
capable of 25 kv. The PCA is usually operated at 11 to 14 kv with deuterium at 8 
torr. At 13 kv the energy storage is 3.5 kj , the peak current 500 ka. The half 
cycle time is 2 y sec. 

Figs. 1, 2, and 3 are all taken at a time when the discharge current is flow- 
ing from the face of the center conductor to the outside conductor. 
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Hydromaqnet i c Heating of Plasma 
For Production of Controlled Fusion Reactions 

Robert A. Gross 
Columbia University 
New York, N.Y. 10027 

This research is directed at creating a high temperature, 
fusion reacting gas by means of a very strong shock wave, and 
understanding the resulting plasma properties and behavior. The 
shock speed required to create a thermonuclear reacting plasma 

Q 

in room temperature deuterium gas i s 1 to 3x10 cm/sec, with 
corresponding plasma temperature of 1 to 10 keV . 

A small, one meter long, coaxial, electromagnetic shock tube 
has been used for many years at Columbia University to study 
ionizing shock waves in hydrogen at speeds up to about 2x10^ cm/sec 
(T~40 eV). During the past year, this device has been used prin- 
cipally for developing better plasma shock diagnostic techniques, 
and, in particular, laser phasographic and holographic nanosecond 
photographic methods. 

A higher energy electromagnetic shock tube was completed in 
1969 .- This new device is 3 meters long and has an outer diameter 
of 22.5 cm. Its shock waves are driven by a 2x10° amp. current 
delivered from a 3x10^ joule, 100 keV capacitor bank. A preshock 
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magnetic field up to 10 kG . is applied usually in the azimuthal 
plane. This high energy shock tube is now operating and some 
preliminary results obtained in it are: 

1. The experimental results are very reproducible. 

2. Ionizing shock speeds of 70 cm/psec have been created 
in hydrogen with a preshock pressure of 25 mTorr. 

This shock speed correspondes to an acoustic Mach num- 
ber ~ 550 and an Alfven Mach number ~ 4. A laser inter- 
ferometer measurement showed a shock compression ratio 
of ~ 4 indicating a post shock plasma density of 

n ~ 6x10*^ cm ^ . 

3. At 100 mTorr initial pressure, shock speeds of 30-50 cm/psec 
have been obtained. At 30 cm/psec, the magnetic shock 
structure thickness is 35 cm. The shock wave is followed 

by a current free, hot, uniform plasma about 45 cm long 

(at the 150 cm tube position). This plasma is followed 

by a diffusing contact surface current. The volume of 

4 3 

shock heated plasma is approximately 10 cm . 

4. The shock heated plasma produces copious X-rays. Pre- 
liminary X-ray spectra indicate that at M ~ 400, the elec- 
tron temperature T^ « 500 eV and an energy analyzer indi- 
cates an ion temperature Tj « 300 eV . 

An improved vacuum system has just been installed so that the 

device can be operated with room temperature hydrogen at pressures 

13 - 3 

as low as 1 mTorr (n = 3x10 cm ) with an impurity level less 
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than one percent. Correspondingly higher shock speeds and higher 
plasma temperatures should be obtainable. 

These strong shock wave experiments are being simulated on 
a computer by solving the corresponding shock tube hydromagnet i c 
initial value problem. Comparison between the test data and com- 
puter predictions has provided considerable insight into shock 
formation time, diffusing current layers and other nonsteady phe- 
nomena . 

A theoretical study is in progress of the dynamics and prop- 
erties of the plasma created by a giant laser pulse heating a 
solid deuter i urn- tr i t i urn sphere. Preliminary results indicate that 
for very short, high power laser pulses, the fusion energy released 
exceeds the laser light input energy for laser pulses of greater 
than ~ 7x10^^ ergs. The possible benefits of surrounding the fuel 
with a heavy Z material are being investigated. 

A theoretical study is in progress to determine the behavior 
of a D-T fuel pellet injected into a fusion temperature plasma. 
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Plasma Turbulence Research 



by 

Norman C. Jen 

The City University of New York 

In recent years, one of the most important problems in the 
field of plasma research and the development of thermonuclear 
fusion control is the understanding of the characteristics of 
plasma turbulence. Due to the lack of physical pictures of the 
turbulence and the inadequacy of the non-linear mathematical method 
for describing the complicated dynamics, significant progress in 
plasma turbulence research becomes nearly impossible. The objec- 
tives of this research are (a) to obtain the actual physical 
pictures of the turbulent plasma through experimentation, (b) to 
understand the interactions and energy transmission mechanisms 
among particles and waves by the quasi-linear theory, and (c) to 
develop a computer program for such a problem with the feeding of 
the experimental results as input data. 

A set of systematic experiments has been conducted in the 
specially constructed quiescent plasma device with a magnetic 
field, ranging from 50 to 6,000 gauss. The highly ionized plasma 
column, which is 0.75 inch in diameter by 20 inches in length, has 
a density of 10^*^ - 10^^ particles per cu. cm. at a temperature 
of 2,400 Kelvin degrees. The dynamics of the turbulent plasma has 
been successfully recorded by the newly developed high speed plasma 
camera at one end of the plasma column. The results obtained are 
as follows: 
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(1) At a low magnetic field of approximately 200 gauss, the insta- 
bility sets in and gives a symmetrical pattern of convective cells. 

The convective cell has the characteristics of zero frequency of 
oscillation, and maximum efficiency of confinement. The convective 
cell is shown in Figure 1, 

(2) At a magnetic field of moderate strength of approximately 1000 
gauss, the frequency spectrum shows coherent waves with a predomi- 
nant frequency in the order of 10, 000 cycles per second. The value 
of the frequency agrees with theoretical predictions. 

(3) At a higher magnetic field above 3000 gauss, the plasma column 
with a longitudinal current becomes unstable in the form of 
Kadomtsev helix. At the same time, the plasma column gives a con- 
tinuous frequency spectrum similar to the spectrum in ordinary hydro- 
dynamic turbulence. The photographic pictures show the rotation of 
the plasma column and the pitch of the helix as a function of the 
magnetic field. In addition, the pictures also indicate the con- 
tinuous change of density distribution and the collective oscillations 
of the plasma across the plasma column. This oscillation is the 
primary factor for the continuous spectrum and the most important 
element for the plasma turbulence. The pictures are shown in Figure 2. 

The theoretical studies on the quasi-linear theory lead to 
following results: 

(1) The collision term in the Boltzmann equation, in the form of 
Krook model or Fokker-Plank model, has significant meaning in the 
regime of medium density plasma which is the case in this research. 

(2) The disturbed distribution function of plasma column can be 



66 



3 



excited into a turbulent state or can be transformed into a Max- 
wellian state throuqh the interaction or the energy transmission 
mechanism among particles and waves* 

For a computer program, the time dependent magnetohydro- 
dynamic equations in the Lagrangian system with two space coori- 
nates have been written and successfully tested. This program has 
a number of defense applications. 

By using the newly developed plasma camera technique in the 
quiescent device, combined with the theoretical analyses and the 
Lagrangian computer program, the research on plasma turbulence has 
been productive in many areas and a much better understanding of 
plasma turbulence is likely to be attained in the near future. 




Figure 1 




Figure 2 (B = 5000 gauss) 
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Heating and Accelerating Mechanisms 

In the Plasma Focussed Discharge 

0. K. Mawardlj Plasma Research Laboratory 
Case Western Reserve University, Cleveland, Ohio 



Introduction 



The production of hot dense plasmas by means of focussed discharges has 
known for some tlme.^ The early accounts on plasma focussing investigations 
reported In the literature referred to two distinct electrode configurations 
coaxial accelerators^ and the metal-walled discharges. ^ 



been 

the 



The usual explanation for the formation of a dense focus by a coaxial accel- 
erators has rested on the assumption that the dynamics of the plasma consists of 
an Initial acceleration phase of the plasma In the "barrel" of the gun, followed 
by a fast z-plnch effect. Experimental observation of the plasma focus by means 
of optical framing pictures and x-ray photographs^ seem to Indicate that the above 
proposed explanation for the plasma focus is perhaps correct since It describes 
the appearance of the discharge in its broad lines. Upon close scrutlnlty, how- 
ever, this model either breaks down or is unable to account for several puzzling 
features of the plasma focus. 



The time associated with a hydromagnetic instability is often associated 
with the ratio of pinch diameter to the local Ion speed. When this criterion is 
used it turns out that this time Is almost two orders of magnitude shorter than 
the lifetime of the plasm focus. Even when due consideration Is given to 
shock heating and coirpresslon it Is evident that the simplified picture of the 
linear pinch is unable to account for the observed life time and rate of heating 
of the plasm. 



Other features of the plasm focus which are still not understood are: 

1) The remrkable flatness of the front of the discharge 

2) The mechanism for the neutron production 

3) The shape of the x-ray spectrum 

The sensitive dependence of the discharge characteristics 
on the shape of the electrodes. 

Proposed Research Program 

1) Theoretical Work 

Two problems cire proposed to be Investigated 
a) Eynamlcs of Focussed Discharge 

There are two features of the focussed discharges which are very important to 
clarify. In many experiments the discharge tak:es place in a static neutral gas 
background. This is in contradistinction to the pulsed gas discharges where the 
background gas pressure is very low. Now, it was pointed out by Kuhkel and Gross' 
that the dynamics of an ionizing front resembles that of a detonation wave. This 
is because, the Lorentz forces are sensed by charged particles only. Consequently, 
in an ionizing front the charged particles can be visualized as sensing the electro 
mgnetlc forces as soon as the front overtakes the particles. A parallel situation 
occurs in a detonation front where the energy is released at the front, l.e., as 
the unbumt gas crosses the front it bums thus releasing the energy. An immediate 
consequence, of this analogy is that the ionizing fronts have a limiting velocity ■ 
(the analog of the Chapman- Jouguet criteria). The collapse of the discharge will ^ 
therefore be set by the hydromgnetic counterpart of the Chapman- Jouguet condition.' 
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Another very Interesting feature of the focussed discharge was discussed 
by Morozov. He suggested that the high tenperatures and densities found in the 
plasma focus may be explained on the basis of hydroiagnetic arguments when the 
conponent (B-VB) of the magnetic pressure is retained in the calculations. 

Morozov was able to demonstrate that for purely azimuthal magnetic fields as 
encountered in focussed discharges, conservation of momentum of the plasma 
requires it to be compressed mainly by the (B*VB) terms as it flows towards the 
focus . 1 

Density conpressions of the order of ^ obtained, where M^ 

is the Alfvenic Mach number for the plasma flow. For an initially cold plasma 
(of the order of 1 ev. ) M. can be quite large ('^-10'^). It is readily estimated 
that for this case a proton tenperature of 5 kev can be reached. 

A model for the dynamics of the plasma flow which incorporates both of the 
aforementioned features and is able to account for most of the observed charac- 
teristics viz. , 

a) the formtlon of the bubble, 

b) the flatness of the front of the discharge, 

c) the jet-like character of the focussed discharge and 

d) the sensitive dependence of the appearance of the 
focussed discharge on the shape of the central 
electrode (hollow or solid) 

is the hydromagnetic analogue of the collapsing cavity produced by the penetration 
of a jet in a liquid . 

There are more subtle reasons for the parallel between the two phenomena. The 
discharge proceeding in a neutral gas, the front will becone an ionizing front and 
these fronts are known to be fairly stable and behaving like a detonation wave. 

The use of a hollow electrode, allows the forward moving jet to appear. Actial 
time integrated photograph for this phenomenon clearly demonstrates this effect. 

It is thus proposed to develop the hydromagnetic theory for this model taking 
into account the corrpresslbillty of the plasma. The properties which are proposed 
to investigate in order are 1) the flow pattern, ii) stability of the flow, ill) 
expected heating in the "pinched" region. 

b) Theoretical Study of Electron-Heating 

A characteristic feature of the focussed discharge is the ^pearance of 
discontinuity in the current. This discontinuity is acconpanied by neutron-production. 
It is often stated that the resulting electric associated with this discontinuity 
is responsible for the local acceleration of the deutrons. 

It is reasonable to expect that this strong electric field appearing in the 
focus may lead to onset of Instabilities and turbulence. This in turn causes the 
turbulent heating of the electrons. Since the electrons appear to be colder than 
the ions, ion-acoustic instabilities are strongly Landau damped^ and an other 
instability has to be Invoked to explain the anomalous electron heating. 

It is proposed to examine, the turbulent behavior of the electrons under the 
assunptlon of magneto-acoustic instabilities. The electron heating is proposed 
to be estimated from simple linearized turbulence theory. The tenperature thus 
predicted will be compared with those measured experimentally as described in the 
section on experiments. 
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2) Experimental Work 



In addition to the routine monitoring of the focussed discharge by measuring 
the current voltage, probing of the magnetic field and visual observation of the 
discharge by means of an STL Image converter camera, It is proposed to conduct 
the following Investigations 

A) ivfeasurement of Electron Density 

Average values and spatial distribution of the electron density In the discharge 
are proposed to be measured by means of two separate techniques 

a) Spectroscopic Techniques 

The mean value of the electron density in the focus or elsewhere in the discharge 
will be estimated by the following standard techniques 

1) Doppler shift 

2 ) Line broadening 

3) Continuum and discrete line Intensity. 

The measurements will be made with a McPherson spectrometer. Eight channels are 
available in the Instrument for time resolved spectroscopy. 

b) Laser Interferometry 

A laser Interferometer suited to the measurement of transient plasma density 
fluctuations has been constructed^® through the modification of an Ashby-Jephcott 
interferometer.il The Increased sensitivity results from a technique Introduced 
by Hooper and Bekefll2, but unlike their design our interferometer works well on 
one shot experiments. Our interferometer can respond to electron density fluctu- 
ations as short as 15nsec.l3 

B) Neutron Production Mechanism 

The purpose of this part of the experimental program is to clarify the mechanism 
of neutron production. In particular the aspects that will be Investigated will 
be the total neutron production and their angular distribution. The effect of the 
design parameter (condenser voltage, gas density, electrode configuration and dis- 
charge period) on the neutron flux will also be studied. 

The angular distribution of the neutrons emitted from the focus has been used 
in the past as a criterion to corroborate or disprove the validity of the "boiler" 
theory of Filipov. Although we too will check this angular distribution we feel 
that the mechanism of neutron production may be better understood by examining the 
deuteron energy-spectrum. The neutrons will be detected by means of a neutron 
activation counter similar to the one developed by Lanter and Bannerman.^^ 

The deutron energy spectra will be inferred from a time of flight technique. 

A drift tube will be attached to the coaxial accelerator and the arrival of the 
particles will be detected by a solid state detector and counted with a high-speed 
single channel analyzer (Hewlett-Packaird Models 5583 and 5590) . 
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GENERAL METHODS TO INVESTIGATE NONLINEAR 



PLASMA INSTABILITIES 

F. Cap, D. Floriani, F. Herrnegger, N. Ortner 

Institute for Theoretical Physics 
Innsbruck University, A-6o2o Innsbruck, Austria 



A . Thermodynamic_ 2 otential_and_nonlinear_stabilitY_in_magneto- 

hydrodynamics , F. Herrnegger 

The stability of a plasma system in a stationary equilibrium 
depends essentially on nonlinear effects due to the interaction 
between the variables determining the evolution of the system 
or due to the dependence of the phenomenological transport 
coefficients on the instantaneous state of the plasma. Prigogine 
/!/ introduced the concept of a thermodynamic local potential 
into the equations of an electrically nonconducting nonequili- 
briura fluid. In the following this concept of local potential 
has been applied /2/ to the magnetohydrodynamic Benard problem. 
The applicability of this method, to the investigation of the 
stability behavior of dissipative flows is due to the fact that 
convective terms appear in the functional of the local potential 
besides the dissipative terms (viscosity, heat conductivity, 
etc.). This is the feature which makes it possible to construct 
a thermodynamic theory of hydromagnetic stability. 

The Benard problem consists in a horizontal layer of plasma 
in. which a temperature gradient is maintained by heating it from 
below in a constant gravitational field. 

We are concerned with magnetohydrodynamic systems that satis- 
fy the ideal MHD basic equations including heat conductivity 
and viscosity 111 . 

We construct a negative semidefinite_^function $ by multipli- 
cation of the momentum equation by ^ , the equation of heat 

^o^v 8T ^ 

conduction by — ^ yy , the equation of induction by yy,and by 
adding. the resulting equations the left hand side defines, the 
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function 4> : 
^ = - 



f 3Tx 2 

2 ''St'' 



T 



^,^.2 / 2 
T 

o 



] 



( 1 ) 



This expression contains only derivatives with respect to time 
t and represents a negative semidefinite function which vani- 
shes identically in the steady state and is smaller than zero 
otherwise: 4> i 0 (2) 

We integrate the result (1) over the volume V, use the principle 
of local potential 121 of the steady state (characterized by 
bared quantities q(x)) and together with the condition (2) we 
get the result 



^4'dV = <|){q(x) ,q(x,t) } (3) 

V 

where is a very complicate functional of the time-independent 
quantities q (standing for v, etc.) and of the analogous in- 
stantaneous quantities q(x,t). In the steady state the left-hand 
side of (3) vanishes, this means the first variation of (|) is 
zero. The Euler-Lagrangian differential equations resulting from 
this variational principle are the steady-state equations from 
which we have started. Thereby a variational principle containing 
convection and dissipative terms has been given (to treat the 
MHD Benard problem) . The principle is useful for numerical treat- 
ment of stability problems. 



LY§£yD2Yl§_§5ability_theory_in_plasma_ghysics , p 
N. Ortner 

The classical stability theory of Lyapunov is successfully 
used for mechanical systems with a finite number of degrees of 
freedom. This theory describes the stability behavior of solutions 
of ordinary differential equations (or systems) without knowledge 
of the explicite solutions. These informations are obtained from 
functional-theoretical investigations of the Lyapunov function. 

In order to evaluate the conditions for stable domains of the 
solutions, the knowledge of the Lyapunov function and of the form 
of the differential equations is sufficient# ,• The effective 
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determination of the Lyapunov function is the main problem. How- 
ever, the governing differential equations in plasma physics 
are partial differential equations /3,4/. In order to apply the 
Lyapunov theory to stability problems in plasma physics v:e have 
transformed the MHD basic equations in ordinary differential 
equations by the help of similarity transformations /5/. These 
transformed ordinary differential equations are then the starting 
point for stability search using the Lyapunov criteria /4/. This 
way of treating the stability of a system has the advantage that 
no linearization of the problem has been undertaken (for more de- 
tails of our preliminary investigations on this problem see Ref. 
/5/). 

C • A_general_algorithm_for_the_investigat ion_of _glasma_insta^i- 
bilitles, D. Ploriani 

The treatment of the full non-linear equations of any plasma 
theory implies considerable mathematical effort. For this reason 
approximate solutions are searched by simplifying the equations 
or by using iterative procedures. The usual linearization can be 
interpreted as the first step of such an iteration. On analysing 
this linearisation one easily perceives the close relationship 
with Taylor series expansion in time t 76/ . V/e assume that every 
quantity can be expanded in a Taylor series: 

?(r,t) = ?„(?) t^/n!, 

n = o 

? being the vector with components p (r), v (r), v (r), 

“ ^ ^ ^ x^n 

V T (r), B „(r), ... . Introducing this into the MOD 

z ^ n n X j n 

equations we get recurrence formulas for the coefficient func- 
tions '• 

Instead of differential equations one now has equations contai- 
ning only algebraic operations and applications of the del-opera- 
tor V to already known functions; approximations of arbitrary 
order can be obtained in a systematic manner. If one starts from 
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a configuration near equilibrium the recurrence formulas can be 
simplified. It also seems possible to use this method in statisti- 
cal mechanics /?/• "^he most interesting application may be the 
computation of the initial time behavior of a system near equi- 
librium by help of high-speed computers. 
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RESISTIVE DRIFT WAVES DRIVEN BY A DENSITY GRADIENT 



G . Auer 

Institute for Theoretical Physics 
Innsbruck University, A-6o2o Innsbruck, Austria 



Drift waves and drift instabilities nay arise in a plasma 
whenever gradients exist. Since such gradients, e.g., of densi- 
ty, temperature, magnetic field, etc., occur in any plasma de- 
vice and since drift instabilities have a strong influence on 
plasma losses the problem of drift waves is of fundamental im- 
portance for the problem of plasma containment. 

A suitable device for studying oscillations and waves in a 
plasma is the Q-machine. A great number of experiments in these 
quiescent plasmas was dealing with spontaneous oscillations which 
under certain conditions could be identified as density-gradient 
driven drift waves /!/. The usual theory of these drift waves 
/I, 2/ which is based on a two-fluid model assumes isothermal 
behavior of ions and electrons, i.e. infinite thermal conducti- 
vity. We wanted to examine the other extreme, the case of zero 
heat conductivity and thus of adiabatic oscillations and to see 
what changes are brought about by this assumption. 

The analysis is based on the first two moments of the Boltz- 
mann equation, i.e. the equation of continuity and the equation 
of motion for both, ions and electrons. Finite electron inertia, 
ion motion parallel to the magnetic field, and parallel electric 
resistivity are taken into account /3/* Since we consider a low-B 
plasma variations of the magnetic field resulting from plasma 
motion may be neglected. The electric field is derived from a 
potential (electrostatic waves). As the frequency of the oscilla- 
tions we are dealing with is much smaller than the plasma fre- 
quency quasineutrality of the plasma is a good approximation. 

The equilibrium (zero order state) is assumed to be as follows 
the magnetic field is homogeneous and has only a component in 
the z-direction, a density gradient lies in the x-direction, and 
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ions and electrons drift oppositely directed in the y-direction. 
The linearization of the problem is carried out in the well 
known manner, by superposing a slight perturbation onto the 
equilibrium quantities and neglecting terms of second order 
in perturbations. 

It is convenient to assume the zero order density to have 
an exponential shape. In this case the resulting partial diffe- 
rential equations have constant coefficients and we may look 
for solutions of the form exp i(k v + k z - ut). The condition 
that the determinant of a set of homogeneous equations must 
vanish for a non-trivial solution leads to the desired disper- 
sion relation. It is a polynomial of 7th degree in the frequency 
as a function of the wavenumber. The dispersion relation was 
discussed for zero and for large parallel wave number. It was 
solved numerically for the range of small parallel wave numbers. 
The result of numerical computation is shown in the diagram, i 
fig. la and lb, where we have plotted the normalized phase • 
velocity against the parallel wavenumber. 

The- results we have obtained are as follows: The modes which 
are already present in the isothermal case and which are labeled 
1 to 5 in the dispersion curves are only sightly changed by our 
assumption of adiabatic plasma behavior. They were already dis- 
cussed by P.P. Chen 121 . Besides them two new modes occur: Por 
vanishing parallel wavenumber the mode labeled 6 in the diagram 
propagates with the yfold ion drift velocity (where the adia- 
batic exponent y = c^/c^). The electrons oscillate in the x- 
direction, the ions in the y-direction. The mode number 7 tra- 
vels for k =0 with the yfold electron drift velocity. Por ' 
the limit of vanishing electron mass the perturbations of the 
ion velocity, the density and the electric potential vanish. 

Thus this wave is the propagation of a perturbation of the 
electron velocity and electron pressure, respectively. When 
the parallel wavenumber increases from zero the two modes keep 
up their identity over a large range until they turn over into 
acoustic waves propagating with the adiabatic sound velocity, 
of course. As concerns their growth rate the imaginary part of 
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their frequency is so small compared to the real one that they 
are practically stable. 

The investigations will be continued by inclusion of ion 
viscosity and thermal conductivity. 
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Variation of the normalized wave velocity ej L 
with the parallel wave number I (reoiproonl den- 

•ity deoay length smaller than perpendioular wave nuaber) 




Pig, 1b. Variation of the normalized wave velocity CJ/ojL^ L 

with the p€U»allel wave number i (reciprocal density 

deoay lenirth greater than perpendioular wave number) 
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DISSIPATIVE MAGNETOHYDRODYNAMIC INSTABILITIES 



E.H. Jager, P. Unteregger 

Institute for Theoretical Physics 
Innsbruck University, A-6o2o Innsbruck, Austria 



It is well known that finite electrical resistivity and 
finite viscosity may cause and drive important classes of MHD 
instabilities, the so-called "dissipative" instabilities. As 
an understanding of these instabilities is of crucial importance 
for the mastering of plasma containment a considerable deal of 
work has been done on this problem. 

The combined effects of non-zero flow velocity, resistivity, 
viscosity, of a sheared magnetic field on the stability of an 
incompressible plasma sheet being in a state of steady equili- 
brium flow has been investigated by E.H. Jager /!/. Starting 
from the basic MHD equations, the convective equations for 
gravitational force and resistivity, and from the perfect-gas 
equation of state, we consider an infinite horizontal slab of 
a quasi-neutral, incompressible plasma. The equilibrium magne- 
tic field is a sheared field which depends only on y but is 
oriented perpendicular to the y direction which is defined as 
the vertical direction. The equilibrium velocity v^ depends on 
y but is perpendicular to the y direction as well. The gravita- 
tional force per unit mass, g^, which simulates geometrical 
effects has only a y component depending on y, whereas the 
equilibrium pressure p^ is a function of x, y, and z. Viscosity 
is assumed to be constant throughout the slab, and the Boussinesq 
approximation is used. On linearizing the basic equations with 
respect to the perturbed quantities q^(x,y,z,t), with q^^(x,y,z,t) 
= Q(x,y,z,t) - qQ(x,y,z,t), we obtain a system of linear diffe- 
rential equations the coefficients of which either are constant 
or depend on y only. These equations are satisfied by solutions 
of the form q^(x,y,z,t) = q(y) exp[i(wt + lx + nz)J . Substitu- 
tion of this expression into the linearized equations yields 
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a set of ordinary differential equations which after some te- 
dious calculations can be reduced to an ordinary differential 
equation of 6th order for the y component of the perturbed 
velocity v^. The involved coefficients of this equation depend 
on the equilibrium quantities, the frequency w, and the wave 
numbers 1 and n. They can be found in the original paper /!/. 
Several special cases have been discussed. 

If we assume that magnetic shear, resistivity, gravitational 

force, fluid velocity are constant and viscosity vanishes we 

obtain a fourth-order differential equation for v . In the limit 

of vanishing resistivity we find the so-called "Alfven modes" 

2 2 2 

for a moving plasma, with w = w. = (IB + nB ) /p. For non- 

A X Z 

vanishing resistivity there essentially appear two modes one 
of which is damped by resistivity only. 

How consider a viscous fluid of zero resistivity with a uni- 
form velocity field in a magnetic field of constant shear and 
in a constant gravitational field. In this case the dispersion 
relation yields a mode damped by viscosity. 

Next the dispersion relation for a viscous resistive fluid 
shows the combined influence of constant resistivity and vis- 
cosity. 

Finally, preliminary estimates seem to indicate that a 
variable flow velocity changes the non-dissipative as well as 
the dissipative modes but leaves the growth rates of the dissi- 
pative modes unchanged. A more accurate investigation of this 
subject is under way. 

Jager’s research represents a generalization of the work 
done earlier by P. Unteregger /2/. Unlike Jager, Unteregger 
assumed a static equilibrium with v^ = 0. He showed that the 
hydromagnetic modes of an inhomogeneous, incompressible plasma 
are strongly influenced by the dependence of the Alfven frequen- 
cy on the magnetic field. For the case that the field of the 
equilibrium configuration is perpendicular to the wave vector 
of the perturbation he found out that the magnetic field has 
no influence on all possible oscillatory-mode eigenfunctions 
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of the above form. Unteregger finally pointed out that finite 
electrical conductivity and viscosity give rise to new modes 
which do not appear in ideal magnetohydrodynamics. 
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MICROINSTABILITIES IN ARBITRARY PLASMA CONFIGURATIONS 



S. Kuhn 



Institute for Theoretical Physics 



Innsbruck University, A-6o2o Innsbruck, Austria 

Because of the considerable mathematical difficulties ari- 
sing with the use of the Boltzmann or Vlasov equations in ; 
microinstability research the relevant analyses were for a 
long time confined to geometrically simple plasma configura- 
tions (e.g., one-dimensional plane slab). It has, however, 
been shown /I, 2/ by Taylor, Hastie, and Haas that it is possible 
to treat the problem in a much more general manner and to 
account for all geometrical effects from the outset. Now, though 
the stability theory devised by Taylor and Hastie 111 holds for 
geometrically arbitrary plasma configurations it is specified 
to the case of a collisionless plasma undergoing 
low-frequency oscillations with long wavelengths 
in the direction of the magnetic field. Beside these so-called 
flute and drift modes there can, however, arise also high-fre- 
quency modes as has been pointed out, e.g., by Rosenbluth /3/. 
This now is the point where we tried an application of Taylor’s 
and Hastie 's theory to the range of high-frequency 
oscillations , additionally taking into account a 
collision term which might give us an estimate 
of collisional damping. 

Given an equilibrium described by the distribution function 
F(r,v), the electric potential <t>(r), and the magnetic field §(r ) . 
We require the ratio r*j^/L = X to be much less than unity, where 
r^^ is a typical Larmor radius and L is a typical macroscopic 
length of the system. Apart from this assumption the geometry 
of the configuration under consideration is left arbitrary. On 
this equilibrium we now superimpose small perturbations 



o^ the distribution function and of the electric potential, resp. 




These expressions are inserted into the Boltzmann equation ’,jhich 

is then linearized. The collision term is assumed to be of the 

Krook form -v ,-,f. Now it has proved convenient to introduce 
coll 

"natural" coordinates as follows. On each magnetic line of force 
(labeled by two parameters a, 8 such that § = ( Sa/ 3r )x ( 3 8/3r )) 
one can define in any point (determined by some arc length s 
along the line of force) an accompanying coordinate frame con- 
sisting of e^ = S/B and two unit vectors e^ perpendicular 

to e. and to each other. Then the velocity may be expressed in 

^ ^2 ^2 
terms of the variables m = Vj_/2B (magnetic moment), K = v /2 + 

+ e<t’/m (equilibrium particle energy), (}>= tan ^(e^.v/e^.v) 
(azimuthal angle of Vj_ about §), and o = e^.v/|e^.v| (sign of 
v„ with respect to e^). Introducing x = Ar, V = 3/3x, and ex- 
panding the quantities involved in simple power series in X 
the Boltzmann equation reads in the lowest order in X 



3f 



i ie_ f + i ^ f D'S - i § 4- VS .(oqe, + v ,D. )F = ^ 
0)^0 e o o Boo ^l3K -•■l o d(p 

where D' and are complicated operators /H/. 

The solution to this equation is found to be 

2 „ 



f^ = b^ - b^^d-f' ^ + 0(<|) •-<{) )j exp 



I(a^,<J>) - I(a^,<t>')J 



where b^, b^, a^, and I are complicated functions depending 

on known quantities. Prom Poisson's equation we obtain 



Q( 



0)) = )d^r U^(VS^)^/i^1T - t 21e\^d^vf 1=0. 

> L. O O O • J O J 



stability analysis now requires to find out under what conditions 
this equation may have a solution w lying in the "unstable" lower 
u) half-plane. This question may in principle be answered by means 
of a Nyquist plot of Q(w). By a Nyquist plot in our case we under- 
stand the mapping of a closed contour (consisting of a straight 
line just below the real axis and of a half-circle in the lower 
(11 half-plane with its radius tending to infinity) into the com- 
plex Q plane. 

We were able to prove the following properties to hold ge- 
nerally for Q(u)'/: 
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a) Q(o)) has no poles located in the lower half-plane since the 
only poles that possible arise can be shown to lie above or 
(for the collisionless case) on the real axis. 

b) Q(ti)) has no branch points located in the lower lo half plane 
for this function can be shown to be a simple power series 
in ID in the lower (d half-plane. 

Because of properties a) and b) we are able to apply the 
Nyquist theorem which states that there are no unstable modes 
if the Q contour does not encircle the origin. 

c) The infinite half-circle of the to contour maps onto one 
single point Q(«) located on the real Q axis. This statement 
is proved by some more detailed calculations. 

Finally, to complete the formal part of our work, Q has been 
calculated for the collisionless case 1^1 in which it appears 
to be slightly more concise than in the case of collisions. 

Given the expressions for Q(w), our stability problem is in 
principle solved. We have to substitute the specific parameters 
and functions of the system under consideration for their general 
counterparts in these expressions and to calculate the shape of 
the Nyquist plot. Such an investigation will yield sophisticated 
criteria necessary and sufficient for stability. 

Our present work is dedicated to discussing the physics 
standing behind the method under consideration. 
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OVERSTABLE nODES DUE TO RESISTIVITY GRADIENTS 



F. Herrnegger 

Institute for Theoretical Physics 
Innsbruck University, A- 6 o 2 o Innsbruck, Austria 



In a sheet pinch in a plasma of variable electrical resisti- 
vity the following types of instabilities occur: 

a) the tearing instability in which the sheet breaks up under 
the action of the magnetic field into filaments parallel 
to the direction of the current, 

b) the rippling instability in which the locally variable elec- 
trical resistivity changes the current density. 

A coupling occurs between the transverse variation of resisti- 
vity and the steady current associated with the main magnetic 
field. 

The purpose of this paper is to study these problems /!/, 
especially the effect of resistivity gradients perpendicular 
to the main magnetic field on the generation of overstable modes 
in a current-carrying inhomogeneous plasma sheet which is in 
the presence of a magnetic shear field and of a constant gravi- 
tational field g. It will be shown that overstable modes may 
exist. These solutions are connected with the second derivative 
of the mean resistivity (n") i.e. they depend on whether the 
course of the resistivity curve n(x2) in the plasma sheet is 
convex or concave. Furthermore it is shown that the region of 
neutral stability depends essentially on the resistivity gra- 
dient and on n", the second derivative, except in the case of 
homogeneous plasmas. 

We consider a current-carrying plasma sheet supported by a 
gravitational field where the following assumptions are made: 

a) the main magnetic field depends on X2^ §(x2) = (X2) ,0,B^(x2) J 

b) neglecting the thermal conductivity along the magnetic field 
lines, the perturbation in electrical resistivity, n, results 
only from convection: 3 n/ 3 t + v.Vn = 0 (v flow velocity) 
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c) no zero-order diffusion velocity exists and the Boussinesq- 
approximation is assumed to be valid. 

Starting from the linearized MOD basic equations and looking 
for solutions of the form 



f^(x,t) = f’^(x 2 ) expTiCic.x - wt)"] 

we determine the criteria for overstable solutions. This problem 

can be reduced to a selfadjoint differential equation of second 

2 

order with respect to b (b = B” ~ k B 2 , the prime denotes d/dx 2 » 
s = (ic.§)) 

(n^b’)’ + b I'- k^n^ + nn*’ + iwpn + (s^n/ipw)^= 0 (1) 



Multiplying this equation by the complex conjugate of b, inte- 
grating over the plasma volume, and assuming homogeneous boun- 
dary conditions we get from (1) a quadratic equation for w 
which has complex solutions only when the condition 



I b ’ I ^dX2 

~ 2 |. |2 , 
j n I b I dx 2 



/nn" |b I ^dx 2 
/n^ I b I ^dx2 



< 




/n I b I ^dx2 
/n^ |b l^dX2 



( 2 ) 



is satisfied. Therefore oscillatory increasing or decreasing 
modes exist. This condition shows that overstable modes are 
associated with the second derivative of the resistivity and 
consequently depend on whether the curve n(x 2 ) has a positive 
or negative curvature. Therefore we see that, since the left- 
hand-side of (2) is positive semi-definite, overstable modes 
occur where s is not zero. Thus the absence of overstable modes 
found by other authors is due to the fact that making an expan- 
sion of the solution about the point where s = 0 and breaking 
off after the first or second term, results in losing the effect 
of resistivity curvature. Furthermore for s = 0 no overstable 
modes can exist, as can be seen from Eq. (2). When local mean 
quantities are denoted by <> the condition for overstability 
derived from (2) reads 

|b| 

The inequality (2') contains the unknown function b. In order 
to make quantitative conclusions we have to solve (1). An 



:> . k2 - <n::> 



2^ ^<-> 



( 2 ’) 
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approximate solution of (1) which is also valid in the case of 
high resistivity gradients, provided the wavenumber k is suffi- 
ciently high, reads 



X. 



b(x 2 ) = ic^exp^ 

nq L 



ij (1 + |/q df. 



X, 



^2 

+ C2exp|- ij (1 + |) ]^q del 

' ■ L 



(3) 



/ N , 2 0)U ^ s^ , -3/4 d^ ( -1/4, 

with q(Xp) = - k + — + — and e = q'^ — p(q ) 



in i(Dpn 



dx. 



(c^, C 2 constants of integration). This approximate solution 
holds in a region where q / 0 and where the further condition 



« i 



_ 1... (1) 
2 ' dx« ' 



<< 1 



is fulfilled. When the solution (3) is inserted into the derived 
stability criterion (2) or (2’), then we get particular informa- 
tion as to how the course of the resistivity in the sheet may 
be changed in order to reduce the amplitude of the instabilities. 
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RESEARCH ON THE PRODUCTION AND CONTAINMENT OF 
ELECTRON RICH PLASMAS 
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J.D. Daugherty, J.E. Eninger, G.S. Janes andR.H. Levy 
Avco Everett Research Laboratory, Everett, Massachusetts 



I. INTRODUCTION 

An electron rich plasma is defined as a plasma in which the space 
charge density of positive ions is significantly lower than the space charge 
density of electrons. Under these conditions, space charge effects play a 
more important role in controlling plasma behavior than do thermal and 
kinetic effects. The intrinsic stability shown by electron rich plasmas 
under properly controlled conditions is of fundamental importance to plasma 
physics, and has stimulated a great interest in a number of new applications. 
A Landau damping effect due to the strong electric fields in electron rich 
plasmas has been predicted, and offers new possibilities to stabilize flute 
type phenomena in plasmas. 

The production of such plasmas is accomplished in two steps. First, 
a cloud of electrons is injected into a toroidal vacuum chamber by means of 
the so-called inductive charging scheme; this scheme is described in 
Section II. Second, the ion component of the plasma is formed through 
progressive ionization of the neutral background gas in the chamber. This 
ability to control the density of ions and electrons independently allows a 
new and a different experimental approach to stability problems. 

In Section III, we discuss a number of instability mechanisms of 
importance in electron rich plasmas, their effect on plasma production and 
containment, and how they can be avoided. Section IV contains a description 
of the toroidal devices and experimental procedures used in the present 
research. The basic understanding of injection and containment gained with 
our Mark I apparatus is summarized in Section V, where also the progress 
of subsequent Mark II experiments is presented. Section VI contains 
concluding remarks with an outline of important problems that remain to 
be studied in this field. 



II. INDUCTIVE CHARGING AND EQUILIBRIA 

L n 

The inductive charging scheme^’ , which we used to establish the 
initial electron space charge cloud, is based on the adiabatic transport of 
electrons in a time -dependent (increasing) magnetic field. This principle 
is best described by considering the guiding center motion of an individual 



Three novel devices have so far been invented in the course of the present 
research: 1) the solar flare plasma radiation shield^ 2) a heavy ion 

accelerator^ and 3) a source of highly stripped ions."^’ ^ 
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particle under the influence of crossed electric and magnetic fields, as 
shown in Fig. l_for an infinite cylindrical geometry with a uniform, axial 
magnetic field B, It is evident that a steady magnetic field leads to an 
equilibrium situation, since the drift caused by the radial space charge 
electric field E^. is purely azimuthal. When B increases, a radial motion 
due to the induction field is added to the azimuthal drift, thus causing 
the electrons to spiral inward. An electron beam launched at the cylindrical 
wall will, therefore, be accumulated within the cylinder, producing a 
negative space charge cloud and a potential well. The adiabatic approxi- 
mation assumes that Kq « E^., which implies that the total number of 
electrons in a magnetic flux tube is conserved during the process of 
electron injection. Under these conditions, it is possible to establish 
various radial electron density profiles by controlling the time history of 
the injected current. This knob is of major importance in experimental 
studies of electron rich plasma stability, since geometrical factors are 
critical. 

The principles of inductive charging can be extended to toroidal 
geometries. We have demonstrated, both theoretically® and experimentally^, 
the existence of toroidal electron cloud equilibria without the necessity for a 
magnetic field rotational transform. *' 



III. STABILITY THEORY 



In dealing with the macroscopic stability of an electron rich plasma, 
one can distinguish between effects that are related to the electron component 
alone, and effects associated with interactions of the electron component with 
the ion component. 

Regarding the first of these two categories, we have demonstrated^® 
that the high frequency, short wavelength stability of magnetically confined 
electron clouds can primarily be described in terms of a single, dimensionless 
parameter: q = o) = neme/e^B^. Although stability against short 

wavelength perturbations does not exist in a strict sense for any value of q, 
theoretical^® as well as experimental^ evidence shows that the stability 
required for all practical purposes (such as the familiar Penning discharge 
and the applications studied here) can be ensured if q 0. 1." 



A long wavelength instability has also been identified and analyzed 
theoretically. This is the diocotron or slipping stream instability ^ ^ 
which occurs due to the electrostatic interaction between two surface waves 
propagating on each side of a crossed-field electron beam of finite thickness. 
The diocotron instability exhibits a high growth rate even for low values of q, 
but it can be avoided completely by the selection of electron beam and electrode 
configurations, in which the unstable long wavelengths are suppressed. Re- 
entrant beam configurations, like those in cylindrical and toroidal geometries, 
can easily be stabilized, since the beam circumference is an upper limit to 
the possible wavelengths. 



In a recent paper, ^ T. H. Stix has applied these results to C. T. R. conditions. 

The growth rate for this gyro-resonance type instability is~ l/2 q^^ exp (-Z/q), 
which becomes vanishingly small for q ^ 0. 1, There are important 
devices, however, based on the highly unstable conditions prevailing for 
q ~ 1, e.g. , the microwave magnetron. 
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The diocotron waves are negative energy waves. This implies that 
any perturbing influence that removes energy from the wave will cause the 
wave to grow. The flow of currents induced in nearby walls with a finite 
conductivity, such as in the present experiments, could provide enough 

stive wall 
1/2 b 2 and 

could be serious for devices with useful electron densities (n 0 = lO^^-lO^^ cm 
It has been suggested, however, that there should be a Landau type coupling^ ^ 
between the higher mode diocotron waves and the outer electrons in a radially 
decreasing density gradient which would serve to damp the waves at a rate 
sufficient to overcome the resistive wall instability. 



dissipation to drive an electron rich plasma unstable. This resi 
instability^®’ should have a growth rate which varies as n^^/a 



- 3 ). 



In the second instability category, we have predicted^^ and observed^ 
an ion-electron resonance instability due to coupling between the diocotron 
wave and ion oscillations in the potential well of the electrons. This instability 
is only present in devices where ions are accumulated (it is not found in 
Penning discharges, for instance), but it has been shown that it can be avoided 
by keeping the ion density below a critical limit. If the ion is very massive 
(mi/Z large), the ion motion will involve frequencies far below those charac- 
terizing the electron cloud, and the two systems are essentially decoupled. 

But as mi/Z gets smaller (e,g. through progressive stripping of a given 
element) the ion frequencies increase and the coupling becomes stronger, 
implying the possibility of an instability. Unstable conditions can also be 
approached by slowing down the diocotron wave. This is the primary cause 
of an instability that occurs in the present experiments, when the ion space 
charge density produced through ionization of the background neutrals reaches 
a critical value of 10 - 20 % of the electron density, in good agreement with 
the theoretical predictions. 



IV. EXPERIMENTAL DEVICES 

Figure 2 shows the Mark I apparatus, ^ which is constructed of 
aluminum and can be evacuated to « 10“^ Torr. The major radius of this 
torus is 46 cm, the minor radius 10 cm. The magnetic field is induced by 
a capacitor bank discharge through the field coils shown. The rise time of 
the field is SOjUsec, and on this time scale all magnetic flux enters the 
containment chamber through the indicated injection slot, where the field 
lines get "loaded" with electrons emitted from a hot tungsten filament. 
Accumulation of negative space charge then takes place according to the 
principles of inductive charging described in Section II. Subsequent contain- 
ment of the injected electron cloud is accomplished by "crowbarring" the 
field coils, thus causing the magnetic field to stay nearly steady for a length 
of time determined by the L/R time constant of the coils, typically about 
1 msec. This simple way of producing a relatively steady magnetic field 
was quite adequate for the Mark I experiments, since the containment time 
was found to be limited to — 60 jULsec by the relatively poor vacuum conditions. 

The Mark II torus, shown in Fig. 3, is an all-stainless steel, 10“^® 
Torr ultra-high vacuum system designed for long term containment studies. 

In addition to its superior vacuum, the containment chamber in the Mark II 
device is entirely separated from the vacuum envelope, thus permitting a 
great deal of freedom in selecting the size and shape of this chamber. The 
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present insert has a square cross-section, and high injection voltages can be 
applied, since ample space has been reserved for electron gun structures. 
Observe that the electron gun is placed at the outside radius in the Mark II, 
in contrast to its inside location in the Mark I. The choice of location is 
associated with the effects of finite electron temperature on injection, and 
will be discussed below in more detail. 

A method has been developed for producing a magnetic field with 
both the short rise time ( ^ 5 ptsec), required for good injection performance, 
and a long steady state time duration for containment. This ’’ring-down -up*' 
scheme is illustrated in Fig. 4 while Fig. 5 shows the resulting time dependence 
of the magnetic field. 

The primary diagnostic tool is the image current button, shown in 
Fig. 3, which can be used in different modes, both to measure the total 
injected charge and electric field fluctuations associated with waves and 
instabilities. Typical image current button signals are shown in Fig. 5. 



V. EXPERIMENTAL RESULTS 

Since the initial results obtained with the Mark I torus have been 
published elsewhere^, we shall only list the major achievements of this 
work; 



1) Development of semiempirical scaling laws for electron 
injection. It was shown that the limit on total charge derives from a natural 
energy limitation in combination with an anomalous noise or electron heating 
mechanism. 

2) Experimental demonstration of the existence of theoretically 
predicted^ toroidal electron cloud equilibria. 

3) Stable containment of ion densities of up to 20% of the electron 
density, and identification of the ion-electron resonance instability^"^ as the 
limitation on containment time. 



The most important result of the subsequent experiments with the 
Mark II apparatus is the achievement of high plasma densities in combination 
with l<gng containment times. Thus, we have contained an electron density of 
3 X 10° cm"^ for 3 msec, yielding a figure of merit n^T ~ 10^ cm“^ sec 
(compared to « 10^ cm"^ sec for the Mark I). The longest containment time 
observed is 5 msec. Fig. 5 shows a set of oscillograms from these contain- 
ment experiments. The quiescent period preceding the rapid onset of the 
ion instability shown by trace c) in Fig. 5, is convincing evidence of toroidal 
plasma equilibrium. The diocotron frequency of oscillogram d) is a direct 
measure of the average electron density, and proves that no charge has been 
lost through non-adiabatic processes. 



In the Mark I, as well as in early Mark II experiments we used a 
diode type electron gun. A serious drawback was that this gun had to be 
operated space charge limited in order to give acceptable injection results. 
This condition leads to extremely severe requirements on filament emission 
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vvcn for moderate filament voltages, and results in poor injection efficiency. 
Those problcn'is were circumvented with a new triode gun, whose current 
can be coat rolled independently of the filament potential. 

Contrary to the diode, the triode also launches the electrons with a 
relatively low temperature, i.e. , kinetic energy in the cycloidal motion 
perpendicular to the magnetic field, Wj^ , as well as in the motion parallel 
to the field, W n . This reduction of the thermal energy permitted an increase 
in the electron density in the Mark II by a factor of four. The decision to 
locate the electron gun of this device at the outside radius was based on 
criteria for the generation of anomalous noise'^ (which occurs when the magnetic 
field strengths reaches a certain critical value), in addition to the advantages 
of having a longer filament. Although we were aware of the fact that the 
thermal effects would be amplified with this filament location (Wj^ increases due 
to the conservation of magnetic moment, and W j| due to the conservation of 
angular momentum), there was then no reason to believe that the temperatures 
produced by the diode would be large enough to have an inhibiting effect on the 
injection. At the present time, however, it is clear that the most favorable 
conditions for electron injection should be with a triode gun located at the 
inner radius of a torus with such a shape, that the ratio of outer to inner radius 
is as close to one as possible. 

Finally, the long containment times have been achieved primarily as a 
result of the improved injection efficiency; this has reduced the gas desorption 
caused by wasted electrons, and therefore delayed the onset of the ion 
111.- 1, '.biiity. 



VI. CONCLUSIONS 

We have given a review of our research on electron rich plasmas, 
and presented the latest experimental results in a strongly condensed form. 
Although several fundamental wave and instability phenomena have been 
satisfactorily covered, there remain questions that need to be examined 
theoretically and/or experimentally. In particular, the predicted Landau 
type damping should be investigated in more detail. The overall evidence 
suggests that a number of interesting applications should be feasible as a 
consequence of the somewhat unusual stability properties of these plasmas. 

As far as it can be judged at the present time, two important problems 
have to be solved before long term containment can be achieved. One problem 
is wall outgassing that occurs in the Mark II torus in connection with the 
injection process, however, a straightforward solution to this seems to exist 
(involving p rebomb irdm t*nl of the walls with electrons from a separate gun). 

A somewT.at reduced outgassing is required for the burnout condition, that 
snoulci keep the ion density below the critical level of about 0 . 1 n^, in which 
case the ion -electron resonance instability cannot occur. 

The second problem is of a more fundamental nature and has to do 
with the theoretically predicted but not yet observed resistive wall instability. 
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Whether this wave growth effect exists in the present experiments or not, 
can only be determined subsequently to the elimination of the much more 
powerful ion instability. If the resistive wall instability is found to be present, 
one should investigate different methods to control it. A distinct possibility 
is that the suggested Landau type coupling will provide sufficient damping to 
completely cancel the wave growth. This result would be of major significance 
because it demonstrates that one can take advantage of an intrinsic plasma 
phenomenon for the suppression of flute type instabilities. 
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Fig. 1 This figure illustrates the principle of inductive charging. 

When the nnagnetic field is steady the only electric field is 
the radial field due to space charge, and the only drift 
is the azimuthal drift Ej./B. When B is time varying an 
induced azimuthal electric field appears, Eq - rB/Z, at 
radius r. ^ This electric field yields a radial drift speed 
Eq/B = rB/ZB. This drift is inwards when B is increasing 
and vice versa. Since this radial speed will normally be 
small compared to Ej./B, the electron paths will resemble 
tightly wound spirals. 
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Fig. 2 Schematic diagram of the Mark I torus. 
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Fig. 3 Schematic cross-section of the Mark II toroidal system. 



I 



100 



i 



MAGNETIC FIELD CURRENT 



D.C. FIELD COIL 




Fig. 4 Schematic illustration of the "ring-down -up^* method used 
for producing a magnetic field with both a short rise time 
and a long, steady state time duration. The system consists 
of a D, C. field coil inside of which is located an A. C. coil, 
that is pulsed for one half of a cycle in order to reduce the 
D.C. field to zero for a short period of time. Long term 
containment requires a true D.C. source, e.g. a battery 
bank, but it is more convenient to use a capacitor bank for 
containment times ^10 msec. 




500 /J. SEC. 0.1 /I SEC 



Fig. 5 This set of oscilloscope traces shows the result of an experi- 
ment with the Mark II containment apparatus operating in the 
"ring -down-up” mode of magnetic field programming. In 
this case a "crowbarred" capacitor bank was used to create 
a semi D.C. field, shown by traces (a) and (b). Trace (c) 
is a signal from an image current button, shov/ing the onset 
of the ion instability about 4 msec after the electron injection. 
The frequency of the resulting diocotron wave, shown by 
trace (d), is a direct measure of the average space charge 
density; here ^ 1. 5 x 10^ cm“^. 
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FAST PLASMA HEATING RESEARCH 



M. Kristiansen and M. 0. Hagler 

Department of Electrical Engineering 
Texas Tech University 
Lubbock, Texas 79409 

The research under grant AFOSR-69-1757 is divided in four major 
parts: Fast Theta Pinch Plasma heating, CO 2 laser plasma heating 

simulations, theoretical investigations of solid-fusion plasma 
interactions, and exploding wires in external magnetic fields. 

The major experimental effort is concerned with a very fast 
theta pinch (TEF-PEE 1). The system will have a 120kV, 24kJ capa- 
citor bank and uses stripline feed plates and a dielectric switch 
to minimize the inductance. The figure shows a schematic of the 
system and the firing sequence. 

At the present time the system is being tested for high voltage 
breakdown with only one of the four capacitors in place. With one 
capacitor the quarter period of the current (^/4) is 500ns. With 
all four capacitors in place it will be approximately 700ns. The 
plasma is preionized with a 45W, 120MHz rf generator; a slow rising 
bias magnetic fi^ld can also be applied. Preheating of the plasma 
IS achieved by switching a low inductance, low capacitance, uncharged 
capacitor across :he switch. This causes the system to ring and 
forms an approximately 2ev, fully ionized plasma. The main capacitor 
bank is charged in approximately 1 msec with a swinging Marx gener- 
ator. Present investigations are primarily concerned with equipment 
performance, such as current rise times, peak magnetic fields, switch 
firing time and jitter, etc. Future investigations will be concerned 
with studies of shock compression and with end losses and possible 
end stopping mechanisms. Diagnostic tools include, among others, 
holographic interferometry, fast cameras, neutron and x-ray detectors, 
magnetic and electric probes. 

A SOW, CW, CO 2 laser is being used in plasma heating simulations. 
The basic idea is to locate the plasma inside the laser cavity or in 
a 3-mirror, coupled cavity arrangement. Stacks of sodium chloride 
flats are used to simulate the plasma and provide the cavity loss 
mechanism. Theoretical and experimental investigations indicate that 
both these techniques cause enhanced heating as compared to the sim- 
ple single pass heating scheme. The final success of the coupled 
cavity technique appears to depend upon the ability to design a cen- 
ter mirror with very special properties. The experimental investi- 
gations have been made on a CW basis but it is felt that the results 
can be interpreted in terms of pulsed systems. Current theoretical 
calculations are concerned with the optimum mirror configuration and 
with the effects of refraction and self-focusing. 

A theoretical study of solid-fusion plasma interaction is under- 
way. The motivation for the study is, among others, the problem of 
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fusion reactor fuel injection and also the possibility of creating 
a powerful UV source. The theoretical problem is quite complex and 
the models used in the preliminary calculations appear to be much too 
simpleminded. The present effort is concerned with formulating a 
satisfactory and yet tractable model. The model involves the nonlinear 
heat diffusion equation and includes the criteria for establishing 
strong, thermal shocks in the solid as well as the effects of sur- 
face phenomena. Calculations will be made using both a uniform plasma 
slab model and also the case of injecting the solid both parallel and 
perpendicular to the magnetic field into an increasingly dense plasma. 

The exploding wire experiment employs a 20kV, 4.5^, 4nh capa- 
citor and a low inductance switch and coaxial wire holder system. 

The exploding wire can be placed in an axial magnetic field. The 
purpose of the investigation is to study the effect of axial mag- 
netic fields on the pauses (or ’'dwell” times) in the discharge. This 
system is also comparatively fast and no pauses have been observed 
when the wire is brought to explode at the peak current . The current 
interruption when the wire explodes appears to be very fast. The 
axial magnetic field currently used appears to be too weak (few 
thousand Gauss) to cause more than a minor change in the I-V char- 
acteristics of the discharge. 
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ION HEATING AND CONFINEMENT OF NEUTRAL-FREE PLASMAS 



A. J. Lichtenberg and M. A. Lieberman 

Department of Electrical Engineering and Computer Sciences 
Electronics Research Laboratory 
University of California, Berkeley, California 94720 



The heating and confinement of a dense >> 9 hot-ion 

(Tj^ > 1 keV) lithium plasma are the objects of the experimental research. 
Lithium is chosen because; (1) it has a small ion Larmour radius for a 
given magnetic field; (2) it has a very low vapor pressure at room tem- 
perature; (3) its ion charge-exchange cross section is low for most con- 
taminant gases; and (4) source techniques (contact ionization or photo- 
ionization) are available to provide a consistent and controllable initial 
plasma for subsequent heating and confinement experiments. 

Generation and Confinement 

The basic research facility is a pulsed, magnetic mirror compression 
experiment^ (BME-II) with a peak magnetic field of the order of 50 kG and 
a mirror ratio of 1.5. The field rises to its peak value in 0.5 msec and 
decays on a time scale of 10 msec. At present, the initial lithium plasma 
is produced by a contact ionization source and is trapped and heated by 
adiabatic magnetic compression. 

The source produces a two inch diameter column of lithium plasma of 
density n = 10^ cm”3 and temperature Tj^ = 0.2 eV. Langmuir probes and a 
recently developed microwave cavity technique have been used to measure 
the density and temperature. To guide the lithium plasma into the high- 
field compression region, a large volume, d.c. solenoid, producing a peak 
magnetic field of 800 gauss, has been constructed. 

Operating with adiabatic magnetic compression, radial probe measure- 
ments show that after an initial stable period of 100-200 microseconds, 
radial loss occurs in several bursts and continues for as long as a milli- 
second. Ion energies of 40 eV, consistent with adiabatic compression 
theory, have been measured^ . Work has been completed on the task of 
formulating a set of coupled, first order differential equations describ- 
ing the dynamics of the lithium plasma in the magnetic mirror compression 
experiment . 

The production of a lithium plasma by photoionization is now under 
investigation. An experiment to measure the photoionization, using a 
Xenon flashtube uv source, is now under construction. 

Plasma Heating 

A broad range of heating mechanisms are being studied for the heating 
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of both ions and electrons in the lithium plasma magnetic compression 
experiment. Because of the low initial ion temperature, in addition to 
the adiabatic heating from magnetic compression, another ion heating 
mechanism must be employed to obtain ions in excess of 1 keV. Separate 
electron heating is necessary if the electrons are to be independently 
confined, and may also be required to stabilize the confined ions. The 
heating experiments are also viewed as a means of testing heating mecha- 
nisms in an environment in which there is little background plasma. 

Theoretical and experimental investigations of electron-cyclotron 
resonance heating in magnetic mirrors is a continuing program. Analyti- 
cal and numerical work^j^ indicates that an electron moving in and out 
of resonance due to its longitudinal motion is stochastically heated by 
a single frequency r.f. source. A simple, stochastic theory has been 
confirmed both by exact orbit computations and by experimental results^. 

The experiments have demonstrated that a short pulse (0.2 y sec) of 3 cm 
wavelength microwaves of high power (250 kw) can heat electrons above 
10 keV. Subsequent adiabatic heating creates an electron plasma with 
temperature above 100 keV. A second mechanism, which heats both ions and 
electrons is ’’turbulent" heating which is initiated by accelerating the 
electrons through the ions by means of a ringing discharge. The subse- 
quent electron-ion instability has been found to strongly heat the elec- 
trons and probably heats the ions also^. However, the experiment was. 
performed with a pulsed source, producing a large gas load, causing any 
hot ions to be lost by charge exchange. 

Current theoretical work centers on an understanding of the basic 
resonant mechanism by which adiabatic invariance is destroyed, and the 
mechanism which generates stochastic behavior. The Fermi acceleration 
process is being studied as a simple example exhibiting the transition 
to stochastic behavior. The transition from a purely stochastic region 
to one in which adiabatic regions also exist can be predicted analytically. 
The results of this investigation are being applied to both ion and elec- 
tron heating. An apparatus to produce a modulated high power (10 kV, 2amp) 
electron beam has been constructed to resonantly heat the ions in the 
lithium plasma experiment. The beam serves as a mechanism to allow the 
accelerating fields to penetrate the plasma, in contrast to external 
fields which are subject to electron shielding. The collective behavior 
of the plasma is taken into account in choosing the heating frequency. 

The turbulent heating experiment will also be repeated in the neutral 
free, lithium plasma environment. 
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Kinetic Theory of Unstable Plasmas (Grant No. 69-1697) 

J.E. McCune and J.D. Callen • 

Department of Aeronautics and Astronautics, M. I.T. 

j ' ' j . 

I . Program and Objectives 

The general objective of this study is to provide more extensive 
information than is yet available concerning properties of plasma instabilities 
of the type for which there is as yet no clear-cut means of stabilization in 
devices aimed at plasma containment in the fusion regime. Outstanding examples 
include a variety of ion-cyclotron modes (of special concern in mirror machines) 
and ’’drift-universal ” instabilities (of special concern in toroidal devices). 
Common to these modes is the necessity of a kinetic description (as opposed to 
a fluid or MHD description) even in the linear regime. 

Our current effort is divided into four parts: 1) Examination of the 

integral equation for electrostatic perturbations in axisymmetric mirror 
machines^, with particular emphasis on the effects of finite geometry on 
electron Landau damping^; 2) Study of the linear theory of drift-universal 
modes (not stabilized^ by min-B), with emphasis on determining the complete 
unstable spectrum^ — these results provide the necessary information for use 
in the non-linear kinetic theory of these modes^ ; 3) Evaluation of the 
importance of coupling between electrostatic and electromagnetic modes^ ; and 
4) Analysis of wave-particle energy transfer with emphasis on wave- amplitude 
limiting effects for ion-cyclotron instabilities . Much of this work is in 
collaboration with other groups and is partially supported elsewhere. Because 
of space limitations we describe below some specific results only for item 2) . 

II . Some Unfamiliar Properties of Drift-Universal Modes 

Using the "slab" model of plasma equilibria with gradients perpendicular 
to B , we have obtained specific results^ for the marginal-stability properties 
of collisionless drift modes (driven by density gradient only). No magnetic 
shear or curvature effects are included in these preliminary results. Typical 
results are shown in Figs. 1 and 2 for a proton plasma with T^/T^ = 1 . These 
results extend the earlier results of Krall and Rosenbluth® by providing 
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marginal stability curves over a wide range of 5 application to the 

Wendelstein (WII) experiment is discussed elsewhere^. 

The marginally stable modes divide into three groups (Fig. 2). For small 
k-^a^^ (<0.4) there is a ’’low-density” branch and a ”high-density branch” 

(i.e., increasing the density eventually stabilizes modes of this type). For 
k^af > 0.4, there is only one marginally stable branch. For the latter modes 
the- threshold frequency is drastically modified by departure from quasi-neutrality 
and is roughly proportional to (Fig. 1). Since in this case (ja/k^ (threshold) 

is of the order of or less than the ion thermal velocity, , the usual 

approximation relating drift waves (generically ) to ion-acoustic modes is not 
useful. An accurate analytic formula for the threshold frequency is shown on 
Fig. 1. 

The vastly different properties of the three groups of marginally stable 
drift modes shown in Fig. 2 maintain themselves in the unstable regime 
[k^/e* < fixed k^X^] . Figs. 3 and 4 show the growth rates (y) 

and real frequencies (d^) for unstable modes with ~ 5 and ~ 0.1, 

respectively. On each figure the effect of density (hX^^^) is illustrated. 

For = 5 (Fig. 3) the real part of the frequency takes on its ’’classical” 

value ^ (modified for finite kX ) only at one value of k„ ; this formula is 

Di ^ 

therefore not useful in describing any significant part of the unstable spectrum. 

On the other hand, the estimate of the maximum growth rate (and the k^ at 
which it occurs), given in Ref. 8, is fairly accurate. We also show on Fig. 3 
the behavior of and y for very small k^ . Analytic formulas for this 

behavior are shown on the figure. 

For k^a^ = 0.1 (Fig. 4), by contrast, the real part of the frequency 
and the growth rate take on their ’’classical” values® over a fairly significant 
range of unstable k^ * s , particularly at ’’low” densities (curve (1)) and 
’’moderate” densities (curve (2)). However, no single formula, describing either 
^^^z^ or y(k 2 ) over the full unstable range is yet available. 

The behavior of a k^a^^ = 0.1 mode at very high density (curve (3), 

Fig. 4) is particularly interesting. At the critical (maximum) k^a^^ for 
onset, is very close to (but slightly less than) . As we move into 

the unstable regime drops slightly as y approaches its first maximum. 

Then, at about the point at which (ja/k^ot^ approaches unity, oij- exceeds (jl)^ , 
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elec trons get out of resonance with the wave, and y accordingly drops 
dramatically. The growth rate then recovers, rising sharply to a second ' 
maximum near 0.3o). . This second maximum is almost 30 times the first! This 
is a ”nonresonant” instability, first noted by Kadomtsev and Timofeev^^, and, 
if it can occur, could lead to very large cross-field diffusion. This mode 
requires, however, a very long wave length along B (Fig. 4) and probably can 
occur, if at all, only in certain toroidal devices (for example, in a high- 
symmetry stellarator with a small, but rational, rotational transform). 

More extensive results concerning the properties of drift-universal modes 
will be reported elsewhere. Effects of departure from slab geometry are under 
investigation. 
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Rotating R-Field Plasma for Thermonuclear Power Reactors 
M. L. Pool - The Ohio State University 

An inspiriting approach to a realistic thermonuclear reactor is now 
evinced through experimental results and through computer analyses. 

The former is secured by an auxiliary rotating electric field and a 
revolving magnetic field, both at the ion cyclotron frequency, impressed 
upon the main magnetic field of the magnetic bottle employed. The 
latter is secured by submitting various pertinent boundary and operating 
conditions to a 360 IBM computer. 

The auxiliary rotating electric field is normal to the direction of the 
central axis of the magnetic bottle and is obtained from neighboring 
electrical currents. The auxiliary revolving magnetic field is parallel 
to B^. As an example, the electric field has a value of 10 volts cm“^ 
on the central axis and increases to twice that value 10 cm off the axis. 

The auxiliary magnetic field has zero strength on the central axis and 
increases with distance from the central axis. The radial variation of 
B^ is utilized to put an upper bound on the energy that an ion can achieve 
during its cyclotron excitation life time. 

In the above fields a sample numerical calculation was made of the 
motion of a collisionless deuteron that started initially with 4 eV of 
energy and 4 cm from the central mirror axis. The rotating and revolving 
fields simultaneously heated the ion to a maximum energy of 57keV and 
moved the guiding center of the ion orbit to within 0. 5 cm of the mirror 
axis. From then on the ion's energy oscillated between 7keV and 57 keV 
while the guiding center remained within 0. 5 cm of the axis. 

Motion in three dimensions has been also numerically analyzed. 
Although the magnetic mirrors tend to put the ion ahead in phase relative 
to the phase of the cyclotron frequency of the oscillator, the electric 
field impressed by the oscillator continually tends to rephase the ion in 
its orbit so that the gain in phase caused by the mirrors is nullified. 

The effects of scattering on the motion of an ion have been numerically 
analyzed at each 0. 15 degrees of arc of the ion's orbit. Many small angle 
scatterings contribute to give a single large R. M. S. scattering angle. 

The rephasing action of the electric field during the mean free time 
reduces the net large R. M. S. scattering angle. Consequently, nearly all 
the ions which initially have energies of even a few tenth of an electron 
volt, and therefore have many scattering experiences, reach energies in 
the keV range without being scattered out of the system. In most cases 
calculations were made where the small angles varied between 10 and 15 
degrees of arc and where a scattering event occurred whenever the 
deflection probability reached or exceeded 40%. 
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Experimental results have been achieved with a gaseous ion source 
and with a disk plasma focus. Experience with an ion source initiated by a 
laser is currently very limited. However, data tend to point out many 
accrued advantages if a laser source were used in the present equipment. 

In summary the salient characteristics of the current and / or projected 
ensemble are: 

1. Rotating electric accelerating and revolving magnetic guiding 
fields are furnished by a two phase oscillator operating at the 
ion cyclotron frequency. 

2. Ions are given arbitrarily chosen maximum energy limits 
above which no further energy is added. 

3. The energetic ions are maintained in a parking orbit. 

4. Scattering losses in energy, in phase, and in direction of motion 
are essentially nullified. 

5. A laser initiated ion source would provide a copious supply of 
low energy ions for subsequent acceleration in an exceptionally 
clean system. 
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BROKEN SYMMETRIES AND CTR CONFINEMENT 



Daniel R. Wells 
University of Miami 

In 1960 the author began investigations of plasma produced by a modified 
theta-pinch coil^. The coil was contoured into the shape of a cone with a 40° 
aperature angle. He observed that plasma blobs were formed which appeared to 
have a well -organized current and flow structure. These structures were fired 
down a long solenoidal drift tube and their properties were studied in some 
detail^. It was observed that they carried trapped azimuthal as well as longi- 
tudinal magnetic fields. Further investigation revealed that they were extremely 
stable and resisted efforts to disrupt their structure. The densities in these 
plasmoids were of the order of 10^^ particles per cubic centimeter and the tempera- 
tures involved were of the order of 100 electron volts or less. The author put 
forth a heuristic theory describing these structures as plasma vortex rings or 
spheroids and began to make educated guesses about the details of their geometries. 
In the meantime, other investigators became enthused about the concept of plasma 
vorticity and began to look for similar structures in other plasma configurations 
and plasma guns. In particular, Dr. Winston Bostick reviewed his work on plasmoids 
which had been in progress since the early 1950's and came to the conclusion that 
the structures he had been studying were also plasma vortex structures, both of the 
type observed by the author and of other types never before observed or studied^. 

The author decided that since the structures were very stable that it might 
be interesting to attempt to build a thermonuclear reactor based on the principle 
of the interaction of plasma vortex structures. His studies indicated that the 
most useful arrangement might be the injection of a set of vortex rings produced 
in the throats of a set of mirror coils and allowed to interact at the center of 
the mirror geometry‘s. The idea was that since each ring carried a trapped magnetic 
field, in the same direction as its partner approaching from the other end of the 
machine, the resulting trapped fields would reverse the field of the main mirror 
coils and produce a minimum § geometry. The resulting minimum ? confinement would 
las^t as long as the vortex rings continued to exist. Thus the rings would have 
the dual purpose of producing the confinement geometry and heating the plasma. 

The heating effect would be produced by the mutual annihilation of the vortex rings. 
This was the original application of the concept of "vortex burning". 
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A machine based on these principles was built at the Plasma Physics Labora- 
tory of Princeton University in 1963. The machine performed exactly as predicted. 
The field was reversed, the plasma was heated at the center of the mirror and the 
resulting overall structure appeared to remain intact for several hundred micro- 
seconds. High density plasma remained in the mirror for times of the order of 
milliseconds, although this afterglow plasma appeared to be cold. The experiment 
was discontinued at Princetonin 1964 because no one except the author and Dr. 
Winston Bostick at that time believed that the concept of plasma vorticity could- 
possibly play an important role in the CTR program. In 1969, many papers were 
published which attempt to show that "convective cells or plasmoids" are the basic 
cause of "pump out". 

The machine was moved to the Plasma Physics Laboratory of the University of 
Miami where it has slowly been reconstructed and improved, as funding allowed, in 
an attempt to at least bring the experiment back to the densities and temperatures 
attained at Princeton in 1964. This has finally been achieved and the resulting 
structures are now under study utilizing the more sophisticated diagnostics made 
possible by a six year improvement in the state of the art. 

In the meantime, other investigators at General Electric Corp.^, Grumman 
Aircraft and Engineering Corp.^, and lately at Sandia Corp.^, have obtained 
similar results and the concept of plasma vorticity is now a respectable idea. 

The work of Cowan at Sandia is especially notable since he has used explosive 
generators to produce the plasma vortex rings. The generators are equivalent to 
a multi -mega joule capacitor bank and the rings carry trapped currents of three 
to four hundred thousand amperes. The same experiment will shortly be repeated 
utilizing the author's geometry of two rings meeting inside a magnetic mirror, 
only now one will have 2 mega joules of energy in each end of the machine. Cowan 
has also introduced the concept of using a planar shock wave passing through a 
theta-pinch coil, parallel to the axis of the coil to produce the vortex rings. 
This allows him to use explosive generators in producing both the planar shock 
and the theta-pinch field and augurs well for a more conventional machine at high 
energies because it eliminates the critical problem of the inherent coupling in- 
efficiencies of a conical theta-pinch gun. It further eliminates the problem of 
impurities produced when the collapsing current sheet leaves the wall of the 
generator. 
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During the six years that have elapsed while attempts were being made to 
bring the machine up to high energies using surplus and outdated equipment, the 
author, along with Dr. Joseph Norwood, Jr., has been building a theory of plasma 
stability which not only describes the plasmoids generated in the conical theta- 
pinches, but also those generated by Dr. Bostick in his coaxial guns®. The 
theory is entirely new and unique in that it is a true global nonlinear theory 
and not just a small extension of the quasi linear or normal mode cascade theory 
now common in the CTR program. From rather heuristic beginnings, the theory has 
grown to one in which the concept of the space-time and various gauge symmetries 
of the magnetohydrodynamic flow field are utilized in order to find the constants 
Qf the motion involved®. For every symmetry there is a corresponding conserved 
integral (Noether's theorem). If t+iat symmetry is broken, then that integral 
muH not be applied* to the problem. Once all of the symmetries and constraint 
integrals have been found, then one applies the principle of least constraint. 

This principle states that if the total energy of the plasmoid is varied, sub- 
ject to a set of constraint integrals on the flow, then the fewer the number of 
constraint integrals applied, the more stable the resulting structure. Thus, i 
the more symmetries that are broken, the more stable the corresponding structure. 
The minimum number of constraint integrals which apply to a linear (superposable) 
flow field results in a set of flow equations which describe a force-free collinear 
vortex ring. If linearity is sacrificed, then fewer constraints can be used and 
the resulting structure is a quasi-rigid rotator centered on the magnetic guide- 
field lines. Many other types of structures are possible depending on what 
symmetries are broken and thus what constraint integrals must be applied. It is 
possible to have several different types of stable vortex structures existing in 
the same plasma machine at the same time. The plasma will always try to form 
these structures if it can. 

It is interesting to note that preliminary studies of Tokamak indicate that 
it consists of a force-free torroidal shell of plasma surrounding a torroidal 
quasi-rigid rotator. 

Future plans are a projected study and correlation of the measured details 
cTf^tPie structures and the spectrum of~plasmoids predicted by the symmetry 
principles. If current experiments on the colliding vortex ring configuration 
further prove this geometry as a stable configuration, then a much more sophisti- 
cated version of the machine should be built. 
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Computer Experiments on Inertial Confinement 

* 

by C. Barnes and D. A. Dunn 

The objective of this project is to understand the behavior of the 
electrostatic-inertial confinement device^ a proposed method of fusion 
plasma generation in which no dc magnetic fields are applied. The device 
consists of a spherical chamber into which many high energy ion beams are 
injected radially inward. Electrons are created at the outer radius with 
essentially zero energy. It is then hoped that a sufficiently dense^ ener- 
getic plasma will foi’m in the spherical chamber for efficient fusion power 
generation. A report of theoretical and experimental work on this problem 
has been given by Hirsch.^ 

The work performed under this contract is based on the use of numeri- 

2 

cal simulation techniques such as those of Dunn and Ho and of Burger.'^ 

The evolution of a simplified model of this device was simulated on a 
digital computer, A one-dimensional model was used, with only radial 
fields and motion allowed. The ions were injected at a single high velocity, 
radially inward. The electrons were injected at the same radius as the 
ions with a negligibly small energy. 

There are two possible steady state solutions to the one dimensional 
spherical problem. If the computer run is started with the ion space charge 
dominating, an intermediate solution forms. Figure 1 illustrates the inter- 
mediate solution. Shown are velocity vs radius plots for electrons and ions 
as well as a potential vs radius plot. In all three plots the sphere center 
is on the right side and the injection radius is on the left side. The 
intermediate solution is characterized by a virtual anode at a non-zero 
radius and a virtual cathode within it at or near the center. As electrons 
are trapped, they tend to neutralize the ion beam and the virtual anode 
moves toward the center. At low ion beam currents, the virtual anode 
reaches the center and one of the two steady state solutions forms. At 
higher ion currents, the high energy ions streaming through the center 
cause a second virtual anode to form there and subsequently move outward 
until it merges with the first virtual anode. Figure 2 shows such a case 
with the second virtual anode having moved part way out to the first virtual 
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anode. The steady state solution thus formed (type-two solution) has a low 
energy ion plasma inside the virtual anode. The other steady-state solution 
(type-one solution) is one in which the ions stream through the center with 
no virtual anode forming (Figure 3) • 

The emphasis of the work performed this year has been to understand 
the above described behavior in terms of the static and small-signal theory 
of simpler systems of beams in finite regions. The basic model studied 
was a planar, one-dimensional diode into which a single velocity beam is 
injected. There is also a background of oppositely charged particles with 
little or no energy which may only partially neutralize the beam, A simple 
theory has been developed for a partially neutralized beam with an infinite 
mass background, and computer simulations of this configuration agree with 

the theory. An extension of the finite mass background theory of Faulkner 

4 

and Ware to backgrounds of smaller mass than that of the beam was made 
and computer simulations were made to check this theory. It was found that 
two types of instability were present. One was the beam instability of 
Birdsall and Bridges with growth at zero frequency. The second instability 
was associated with the background species and predicted growth near the 
background species plasma frequency. The computer simulations confirmed 
the theory and demonstrated the behavior in the large amplitude, non- 
linear regime. 

The results of these planar diode studies were then used to explain 
the behavior of the spherical inertial confinement device. Tlie conditions 
for existence of the Type-One steady state solution were demonstrated and 
the formation of the second virtual anode was explained. 
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BEAM PROBE MEASUREMENTS OF PLASMAS 



R. L. Hickok and F. C. Jobes 

Mobil Research and Development Corporation 
Princeton, New Jersey 08540 

The density, space potential, change density (P), current 

density electric and magnetic fields associated with 

a plasma can all be measured with a heavy-ion beam-probe. They are 

measured by measuring and mapping the number, energy, and momentum 

2 + 

of doubly charged ions (such as Cs ) created by collisional ioni- 
zation of a beam of singly ionized heavy ions (Cs"*^). Well resolved 
beam-probe maps of an arc plasma show that n, o, and J^ are all 

located in specific and somewhat different locations in the plasma. 

2 + 

The secondary (Cs ) ions are able to provide a localized 
measurement of the plasma because, at the point of creation, they 
undergo a step change in radius of curvature in the magnetic field, 
and so are separated from the main beam and also from ions created 
at other points along the beam path. A suitably placed detector 
thus can monitor the plasma at any arbitrary point, and the entire 
plasma can be mapped by sweeping the beam and by moving the 
detector; instead of moving the detector, however, an array of 
detectors can be used or the beam energy can be varied. 

The basic measurements made on the secondary ions are: 
the number of 2+ ions relative to the number of 1+ ions; the 
difference in energy between the 2+ and 1+ ions; and the component 
of momentum parallel to the main magnetic field which is transferred 
to the 2+ ions. These three measurements give n (actually n*F(T)), 
i and an approximation to A^, the component of the vector potential 
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parallel to the main magnetic field. The electric and magnetic 
fields, and the charge and current densities are obtained from maps 
of i and by taking appropriate spatial derivatives: gradients 

to get the fields, and Laplacians to get the charge and current 
densities. 

Our heavy-ion beam-probe system consists of a Cs^ ion gun 

capable of emitting a 1 ua 5-30 KeV Cs"*^ ion beam, and a combined 

energy and momentum analyzer. The ion gun consists of a thermionic 

Cs^ emitter, an electrostatic lens system, and electrostatic sweep 

plates. The beam from this gun has a spot size 4 mm diameter at 1 

meter, and an energy spread of less than 1 eV. The energy analyzer 

is a parallel plate electrostatic analyzer with the accelerator 

voltage on the positive plate, and a correction voltage on the 

negative plate. The ions enter the analyzer through a slit in the 

negative plate, are deflected by the field in the analyzer, and 

leave through a second slit to be collected on a split detector. 

The difference between the signals to the two halves of this 

detector is greatly amplified and applied as the correction voltage 

2 + 

on the negative plate, thereby maintaining the beam of Cs ions 

centered on the split detector: the required correction voltage is 

the space potential $ at the point where the 2+ ion was created. 

The analyzer is mounted so that the ions are deflected in the x-y 

planej ion motion along the z direction does not effect the energy 

measurement. The detector, however, is split in the z direction as 

well as in the x direction required for energy analysis, and a pair 

of sweep plates near the entrance slit of the analyzer serves to 

center the ions in the z direction. The voltage on these plates is 

proportional to the momentum P acquired by the 2+ ions from a 
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current flowing in the plasma. The sum of the signals to the four 
quadrants of the detector is proportional to the density at the 
creation point. The analyzer is thus able to simultaneously 
measure n, and A . 

A set of four beam-probe maps (n, and 0) of an arc 

plasma are shown in Figure 1. The space potential and density 

maps, although similar in appearance, are not the same; the peak 

of the space potential (~100 V) lies outside the density peak 
13 -3 

(~3 X 10 cm ). The grey scale in the left-hand charge density 

map was chosen so that an excess of electrons is white, an excess 

of ions black, and any shade of grey is neutral. In the right-hand 

map of p the contrast is chosen to show the variations in P: the 

positive and negative peaks correspond to excess ion (Z=2) and 

8 “*3 

electron densities of 4 x 10 cm . The semicircle of electrons 

has the same radius as the arc cathode. A density and current 

density pair of maps is shown in Figure 2. Density contour lines 

at 30%, 50%, and 7 0% are also shown on the J map. The current is 

z 

quite sharply peaked and is located in a region where n*F(T) is 
quite low. 

At the request of the Princeton Plasma Physics Laboratory 
we are planning to install a heavy ion beam probe system on the 
Model T.S. Tokamak. This device, which is being constructed by 
converting the Model C stellerator to Tokamak geometry, will be the 
first Tokamak in the U. S. and is scheduled to go into operation in 
May, 1970. The heavy ion beam probe, which we hope to have 
operational by June, 197 0, will be one of the major diagnostic 
systems on the T.S. Primary emphasis will be in measuring the 
plasma current distribution and the plasma space potential. ^ 

I 

1,1 

Vi 
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The T.S. Tokamak will have a major radius of 109 cm and a 

minor radius of 14 cm. The beam probe system is being designed to 

operate in the regime where the primary toroidal field on the axis 

of the donut is 2 0 kg and the total plasma current is 40 kAmp- 

13 -3 

Expected plasma density is of the order of 10 cm and electron 
temperature a few hundred eV. Anticipated energy confinment time 
is about 50 msec. 

The heavy ion probe will consist of a 10 a amp, 200 kev 

4* -f" 

beam of either Ti or . Measurements will be made right next to 

the plasma limiter, which is in the pumping port cross. The 

primary beam will enter through the top port of the cross and the 

secondary ions will exit through the side port and the pumping duct. 

The secondary ion detectors will be located approximately 200 cm 

from the center line of the plasma. Orbit calculations show that 

with the planned port facilities, it will be possible to map the 

current density, space potential, and particle density from the 

outer edge of the plasma to just inside the center line. 

The anticipated spatial resolution for an individual 

2 

measurement is about .03 cm . For the current density and particle 
density measurements the spacing between grid points can be made 
less than 1 cm if desired; but, at least initially, it is planned 
to use a somewhat coarser grid network--probably 2 or 3 cm. The 
grid spacing for the space potential measurements is a somewhat 
more difficult task because it appears that it will be necessary to 
use multiple electrostatic analyzers, and we have not yet worked 
out the geometric details of this installation. The temporal 
resolution for individual measurements will be of the order of 1 u 
sec, and it should be possible to make a 100 point map in 1 msec. 
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Fig 1 




Fig 2 
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TUNED LASER DIAGNOSTICS 



Dr. R. M. Measures 
Institute for Aerospace Studies 
University of Toronto 
Toronto, Ontario 
Canada 



ABSTRACT 



Preliminary experimental results have been attained from the experi- 
ment designed to test the concept of selective excitation spectroscopy. A 
low pressure potassium arc was used with a tunable ruby laser. The results 
obtained are encouraging for a signal has been observed which can only be 
accounted for on the basis of intensified spontaneous emission. Moreover the 
duration of this pulse is foiind to depend strongly on the electron density. 

REVIEW AND OBJECTIVE 



The concept of using the process referred to as selective excitation 
spectroscopy to diagnose a plasma was first suggested by the authorv^/. This 
process relies upon the ability of a suitably tuned, powerful laser pulse to 
momentarily redistribute the population between the two states of the atoms 
or ions of interest. Under these circumstances the spontaneous emission 
emanating from the upper level of the pumped transition is considerably in- 
tensified and should be easy to identify against the background radiation. 

The duration and amplitude of this pulse of intensified emission are related 
to the local properties of the plasma. It was shown in reference 1 that for 
the limiting case of step excitation by a high power beam of radiation the 
peak intensification of the spontaneous emission depends only on the ratio 
of the populations in the two levels prior to irradiation. It was also 
shov/n that the decay of this pulse would depend strongly upon the electron 
density whilst being weakly dependent on the electron temperature. These 
results indicate that the local conditions within a plasma might be evaluated 
from relatively sin 5 )le measurements of the intensified emission resulting from 
suitable excitation. 

The objective of the cxarrent programme is to test this concept and 
ascertain the possible range of conditions over which it could be applied. 

A study is being made of the intensification of the spontaneous emission at 
6911 °A in potassiinn resulting from the excitation of the 6939 °A line by a 
pulsed ruby laser. An in^jortant by product of this research has been the 
development of a new concept for preionization in the MHD power generation 
field. This involves enhancing both the electron temperature and number 
density of an alkali -seeded pl.asma by pumping the .resonance lines of the 
alkali atoms with a suitably tuned powerful laser 

APPROACH AND RESULTS 



A more realistic theoretical model has been developed to calculate 



131 



the intensification of the spontaneous emission resulting from the selective 
optical excitation of a specific transition within the atom of interest. The 
teiip)oral variation of the exciting laser radiation was represented by an 
analytical expression which was fovuad to agree quite closely v/ith the shape of 
the experimental laser pulse. Both the close coupling between adjacent energy 
levels and the degree of overlap between the laser v/avelength and the v/avelength 
of the line being excited has been taken into account. 

Figure 1 illustrates how the pulse of intensified spontaneous emission 
at 6911 Ain potassium mi.ght “be expected to vary as the laser is progressively 
detuned from the 6939°A line. The conditions assumed in this calculation are, a 
potassivmi plasma with an electron density of 10^^ cm“3, electron temperature 
of 2500 and a laser beam of lo5 watts cm~^. It can be seen from this figure 
that provided the laser v;avelength falls v/ithin three Dopple widths of the 
line centre the pulse of intensified emission becomes essentially independent 
of the degree of overlap. It is also possible to interpret these curves as 
indicating the variation in the intensified emission for different laser 
powers when the laser v/avelength coincides exactly with the line centre. 

Under these circumstances it is found that for laser pov/ers greater than about 
500 v:atts cm"2 the pulse of intensified emission becomes practically inde- 
pendent of the laser power. In Fig. 2 the variation of the pulse of intensi- 
fied emission at 6911 °A is calculated assuming a laser power of 30 watts cm"^ 
at a wavelength coincident with the centre of the 6939 °A line. The electron 
temperature was taken to be 3 OOO K. 

The experimental prograr/irae uses a specially designed Q- switched ruby 
laser and a low pressui'e arc. The lasing wavelength of the ruby is varied by 
controlling its temperature. This is achieved by circulating cooled nitrogen 
within the laser cavity. Figure 3 presents the tuning curve and as indicated 
in the small insert a split Fabry-Perot ring system is used to determine the 
degree of overlap of the lasing v/avelength v/ith that of the 6939 A line in a 
potassiiua lamp. The discharge chamber had to be carefully designed and tested 
to ensiu'e that laser scattered radiation v/as kept to a minimum. 

t 

The first results of this experiment have been attained. They clearly 
establish that the pulse of intensified emission can be identified against 
the background radiation of the arc and that scattered laser radiation is no 
problem. Several tests v;ere conducted in order to ensure that the observed 
signal, an example of v/hich is illustrated in Fig, ^f(a), vas due to intensi- 
fied spontaneous emission. The radiation at wavelengths betv/een the l.aser 
v/avelength at 6939 A and 69 II A v/as monitored and found to have dropped 
effectively to zero by The laser v/as also detimed by about 0.2 A 

where upon the signal at 69 II A disappeared. To determine the degree of 

overlap bet^/een the laser and the 6939 ^A line, part of the lasers output v/as 
fed to a Fabry-Perot interferometer where its rings v/ere con^ared to those 
created by the 6939^A emission from a potassium lamp. The Fabry-Perot ring 

systera for the signal shov/n in Fig, ^i(a) is illustrated in Fig. 4(b), 

A comparison betv/een the theoretical curves and this experimental 
signal js also shov.-n in Fig. 2, The variation of the decay times for the 
691 l^A signe.l against the electron density as measured by a Langmuir probe 
is shovni in Fig. It is found that these results agree with those pre- 
dicted theoretically if the cross-sections, used in the quenching of the 
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overpopulated 6 are taken to have a value of about (l/ 7 ) of that 

given by the Gryzinsky model. This could account for the extended duration 
of the experimental pulse compared to the theoretical profile since the curves 
shown in Fig. 2 were for the quenching rates given by the Gryzinsky cross 
sections. The large experimental error bars in the electron density, indicated 
in Fig. 5 5 sore due to the uncertainties in the Langmuir probe measurements. 

FUTURE PROGRAMME 



The experimental results achieved to data are very encouraging although 
they are preliminary in natiire. With suitable refinements more definitive 
measurements should be possible in the near future. It is clear that much 
work remains to be done in order that the range of applicability of this new 
diagnostic technique can be ascertained. 
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WDWnON fta ABSORPTION LIfC CENTRE, m THE LIfC IS ASSUCD ID BE pOfFlER BRWBCD 
WIH Hm^INIENSITY m WIDTH, /W THE LASER FOCR, O^x) - i0> #\TTS Cr2, 




FIG. 2. VARIATION OF INIENSIFIED EMISSION PULSE AT fflUOA PDR DIROENT EUCITCN DENSITIES 

eajgoA uft in a potasshu piAsm sEuaivav excited ey a pulse of um radiation 
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An Experimental Study of Electrostatic 
Plasma Confinement* 

J. T. Verdeyen, B. E. Cherrington, T. Dolan, D. Meeker 
University of Illinois, Urbana, Illinois 61801 
The concept of inertially confining electrons and ions as 
originally proposed by P. T. Farnsworth and investigated by R. L. 

Hirsch is the subject of our experimental and theoretical in- 
vestigations. Before launching into a detailed discussion of our 
experimental efforts, it is best to review the elementary theory 
of the device as enunciated by Hirsch. 

One first considers a space charge limited emitter of ions 
of mass M which are attracted inward by a highly permeable grid 
structure. We will limit our discussion to the case of cylindrical 
geometry with injection inward towards r=o. If sufficient ion 
current is injected into the interior of the cylindrical grid, the 
potential is raised to a sufficiently high level to repel any 
additional charge. The virtual electrode thus formed on the interior 
has all of the characteristics (in theory) of the original emitter. 

Now if we allow the permeable grid to also be an electron emitter, 
the charge structure of the interior becomes considerably more com- 
plicated. First, the electrons are attracted inwards by virtue of 
the excess positive space charge. Their inertia can carry them through 
the positive space charge and thus buildup a net negative charge. 

Such a picture is quite consistent with a simple solution to Poisson’s 

* 

Much of the work reported here is part of a Ph.D. thesis to be 
submitted by T. Dolan. 
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equation for two charge species with one dimensional radial motion 
and monenergetic particles. It could also be consistent with the 
experimental results obtained by Hirsch at ITT, If one integrated 
the equations inward from the grid, an infinite number of such 
charge layers or virtual electrodes would be formed. Obviously such 
a picture is much too elementary. 

Last year at the AFOSR contractors meeting in Illinois, we 
heard the results of a computer study of such a problem by C. Barnes, 
Although one dimensional radial motion was still retained, no 
restriction was placed on the distribution function of the particles. 
More important, he was able to follow the time history of the 
development of the potentials. His results seem to indicate that 
one virtual electrode or a potential well could be formed, but 
that the interior of this would be filled with a quasi-neutral low 
energy plasma. 

Our experimental efforts have been devoted to the diagnostics 
of a device which was built to approximate the theoretical model. 

The three diagnostic tools utilized have been: ' 

1. laser heterodyne diagnostics using a He:Ne laser 
operating at 640lS 

2. Collimation studies of the excitation light emitted by 
the plasma. 

3. Studies of electron beam deflection by the potentials. 

Both (1) and (2) have been studied in detail in our cylindrical 
system involving two charge species, whereas we are using (3) to 
study electron injection in spherical geometry. 

In our experimental apparatus, the ions are formed by a pulsed 
discharge between a grid and 6 thermionic cathodes placed symmetrically 
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around the periphery of the cylindrical structure. These ions (or 

electrons) are then accelerated inwards toward a highly open (73-92%) 

grid structure which is 4 cm in diameter. Impact ionization of the 

low pressure neutral background gas (1-lOy of neon) within the inner 

cylindrical structure provides the electron source required for the two 

charge specie model. We initially planned to measure the electron 

density within the device by using laser heterodyne techniques. 

Unfortunately, our densities were below the limit of detection by 

11 -3 

laser techniques (2x10 cm ) so we took advantage, instead, of the 

anomalous dispersion between the laser line in He-Ne at 640lS 

o 

and an absorption line at 6402A from the metastable Is^ state of neon. 
Since neon was the background gas for this experiment, a sufficient 
density of these states were produced by electron impact to allow 
us to measure them with a laser. 

However since the metastable state can only be destroyed by 
collisions either with the wall or with another particle, the 
measurement of the radial distribution of metastables does not 
necessarily reflect the radial distribution of the excitation. 

Computer studies indicate that the initial metastable density 
accurately reflects the spatial profile of the source for early 
times in a pulse discharge, but rapidly smears out due to diffusion. 

The spatial profile of the metastable density and the exci- 
tation light indicates a shallow potential well forming on the 
interior of the center grid. Most significant is the fact that 
the depth of the well appears to be critically dependent upon the 
mesh spacing of the grid wires. Consider two cases with the same 
circulating current and the same transmission factor of the grid 
structure, but with different spacings. The metastable density 
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will be considerably higher on the interior for the smaller mesh 
spacing. We attribute this to the shielding effect of the low 
energy plasma on field lines terminating on the grid wires. 

The detailed theory of the measurements and the results will 
be discussed at the meeting. 
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The energy transfer between the various components of a typical plasma 



system of interest here is illustrated in the diagram: 




Eventually the entire chain of reactions set off by an initial fast ion 
either dies away, or it may combine with fusion reactions to create an even- 
larger population of highly-energetic ions. The possibility of such fusion 
and knock-on reactions linking sufficiently to form runaway chains was 

3 

asserted by Gryzinski some years ago as a possible supernova mechanism. 

His estimates of the slow down rates appear to be low, however, and the 
conditions quoted for supercritical chains to occur do not correspond to those 
obtained from our calculations (as discussed in the next section) even when 
neglecting, as did Gryzinski, the contribution to the stopping power from the 
ion-ion Coulomb interactions. 

Even though the individual reaction chains may die away, the effects on 
the energy balance and energy transfer rate of secondary reactions and knock- 
ons of superthermal particles in a hot electron plasma can still be substantial. 



3 



M. Gryzinski, Phys . Re.v. 115 . 1087 (1959). 
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For example, the alpha from a D-T or D-He reaction (born at an energy of 
about 4 MeV) creates with significant probability fast deuterons of about 
1 MeV energy, because of the large resonance in the alpha-deuterium scattering 
at 2.2 MeV alpha energy. 

The most convenient formulation of the population kinetics problem from 

the point of view of numerical calculations combines a Fokker-Planck des- 

4 

cription of the Coulomb dynamics, as given by Rosenbluth, et al , with nuclear 

collision transport terms to describe the disruptive nuclear collisions and 

5 

fusion events, for which a Fokker-Planck approximation is inappropriate. 

The evolution in time of the distribution functions N^(E,t) in a uniform and 
isotropic slowing down medium is determined by the set of coupled F-P 
equations : 

i . . 2 , , 

1^- (E,t) - (E N^E,t)A"^ (E)) - (D"(E)N"(E,t)) 

at S E s ^^2 

= L. cj N^E,T) + I^(E,t), 

where the c| are effective nuclear elastic energy transport and the and T^ 

have their previous meanings as injection sources and Coulomb-process 
slowing down times for the irt species. The functions D^(E) represent the 

4 

M. Rosenbluth, etal., Phys. Rev. 107, 1 (1957). 

^J. Wyatt, et al., op. cit. 
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dispersive component of the slowing down by Coulomb processes, which is 

generally small for ion energies well above the mean electron energy. In 

the numerical calculations on the slowing down reactions reported here, the 

dispersive component of the Coulomb dynamics was neglected. In deuterium 

and especially in mixed D-T backgrounds at high electron temperatures, or in 

systems in which a major fraction of the background ion distribution has a 

mean energy comparable to the fast particle energies, the dispersive slowing 

terms can be important. Detailed numerical treament of both the ion an 

electron distribution dynamics and evaluations of conventional high energy 

injection energy balances in open geometries, which include these terms but 

6 7 

ignore the nuclear scattering, have been published previously, ' 

The coupling of particle species by nuclear reactions requires that in a 
hot electron plasma with good fast ion containment, all of the five particle 

3 

types: proton, deuteron, triton. He , and alpha be considered simultaneously. 
Calculations with the program OLFIN also have an option to include a sixth 
(neutron) species. This allows one to evaluate what the contributions of 
neutron-produced knock -ons to the population dynamics would be in a 
sufficiently large or sufficiently dense system. 

RESULTS AND CONCLUSIONS 

With no geometrical mirror losses and Coulomb slowing contributed by the 
electrons alone, not only substantial energy multiplication but divergent chain 
reactions were found to be possible, as shown in the Figure for a 50-50 D-T 

g 

T. Fowler and M, Rankin, J. Nucl. Energy, Pt. C, 8, 121 (1966). 

7 

D. Fisher, et al. , Proceedings of Int. Conf . on Plasma Confined in Open- 
Ended Geometries , p. 280 (1967) . 
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mixture. With modest fast ion loss rates corresponding to a mirror ratio of 2, the 



divergent chains disappeared, although substantial energy gains remained at the 
higher electron temperatures. 

Calculations including the contribution of the background ions to the Coulomb 
dynamic friction signficantly alter the results given in Fig. B-4 for electron 
temperatures above 30 keV. Divergent fusion chains are no longer achievable 
without including neutron knock-ons . With these knock-ons , the divergent 
chains reappear at an electron temperature of about 80 keV. Further study of the 
conditions for super-critical chains is being carried out using a more complete 
description of the Coulomb dynamics in the lower range of fast ion energies, 
where the dispersive components of the ion scattering may become important and 
the strong D-T fusion resonance occurs . 
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Relative Energy Release from Fast Triton Slowing 

from 400 keV in 50-50 DT Plasma vs . T . 
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INDUCTIVE PULSED PLASMA CURRENT SHEETS 

C* L. Dailey, TRW Inc., Redondo Beach, California 

C urrent Sheet Diagnostics 

An impulsive current sheet is formed in argon by pulsing current through-, 
an 8" diameter flat, spiral coil. At the pressure of 500 millitorr (1.6 x 10 
atoms per cc.), used for these experiments, the current sheet is diamagnetic 
and less than a centimeter thick. It is formed during the first few tenths 
of a. microsecond while the current in the coil is rising in the first quarter 
cycle of the discharge. No pre-ionization is required. 

The accelerator circuit is shown schematically in Figure 1, and phy- 
sically in Figure 2. The capacitor bank comprises 9 small capacitors having 
a total capacitance of 4 microfarads and storing 288 joules at 12 kv. The 
coil is placed in a glass cylinder having a uniform gas pressure of 500 milli- 
torr of argon. 



The physical processes occurring within the current sheet are investi- 
gated by means of miniature probes that measure electric fields, magnetic 
fields and plasma current density throughout the current sheet. In addition 
a laser light-scattering technique is used to measure the plasma mass distri- 
bution and electron temperature. With these measurements we intend to develop 
a physically plausible and consistent description of the current sheet. 

At the time of writing, the probe measurements are complete. The laser 
measurements are underway, but have not yet produced definite results. 



A representative summary of the probe results is shown in Figure 3. 

Shown here are the plasma current density, the tangential electric field that 
produces the current, the radial magnetic field that is supported by the 
current sheet and the space-charge electric field that accelerates the ions. 
The dashed line in the figure is the electron density inferred from the ex- 
pression 



n = 



-Je " 



e E 



( 1 ) 



where electrons are assumed to be the only current carriers and the force 
due to electron pressure gradient is ignored. 



The curves shown in Figure 3 form a consistent picture of the current 
sheet in which a low density electron current front is followed by a thin 
sheet of high density plasma. However, since the electron density was in- 
ferred from the other data, this interpretation is not necessarily valid. 

In fact, earlier investigation with a V* diameter- inductive accelerator^ 
and with electrode- type impulsive accelerators^^ have shown that when the 
plasma density is determined by some independent method, it is not in agree- 
ment with that obtained from Equation (1). This implies that ions must carry 
some of the current, a situation that seems impossible in the inductive geo- 
metry because the coil is presumably not able to supply the required torque. 

Thus, the laser measurement of density is crucial in order to complete 
the description of the current sheet in the present accelerator and to under 
stand the physics of plasma acceleration in an impulsive current sheet. 
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Thruster Performance Evaluation 



A thruster employing an 8" diameter coil has been built and tested on a 
thrust balance. The primary difference between this device and the accelera- 
tor is the use of a fast-opening valve to discharge gas from a small high- 
pressure plenum to the region in front of the coil. Gas density distri- 
butions, measured by means of a fast ionization gauge, showed that only one- 
third of the gas initially in the plenum was close enough to the coil to be 
accelerated. The efficiency, based on the mass accelerated, was higher than 
that for an MPD ammonia thruster for xenon and krypton propellants and was 
about the same for argon and ammonia. ^ 

For many reasons, e.g., pulsed-mode of operation for spinning satellites, 
possibility of accelerated life testing, system fully operational at atmos- 
pheric pressure for systems integration and launch pad testing, the pulsed 
inductive thruster appears attractive. This initial evaluation of its per- 
formance shows the device to be feasible although development of a better pro- 
pellant injection valve is essential. 

The thruster is shown on the balance in Figure 4 and is seen during firing 
in Figure 5. 

Advanced Thruster Investigations 

The pulsed inductive thruster employs the magnetic pressure between the 
coil and the current sheet to accelerate the propellant mass. In effect this 
field "contains" the plasma away from the spacecraft. 

A different method of accelerating plasma is the "field-annihilation" 
effect. This is illustrated in Figure 6 for a pair of oppositely directed, 
flat, spiral coils. When the axially accelerated current sheets collide 
midway between the coils, the oppositely directed radial magnetic fields 
behind these sheets diffuse into each other and are annihilated, transferring 
their energy to radial plasma motion. For the geometry shown, some plasma is 
driven radially outward and the rest is compressed toward the axis. This 
is convenient for experimental reasons. However, in a thruster, concentric 
coils would be used with plasma in the annular space between the coils being 
ejected axially as the radially moving current sheets diffuse together. 



1. Dailey, C. L. , "Investigation of Plasma Rotation in a Pulsed Inductive 
Accelerator," AIAA J. Vol. 7, No. 1, p 13 (1969). 

2. Lovberg, R. H. , "Acceleration of Plasma by Displacement Currents Resulting 
from Ionization," Proc. 6th Int. Conf. on the Ionization Phenomena in 
Gases (1963). Also General Atomic Rept. GA-4363. 

3. Dailey, C. L. , "Thrust Measurement for a Pulsed Inductive Thruster," 
Special AFOSR Rept. Contract No. F44620-68-C-0042 (1970). 
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PULSED INDUCTIVE ACCELERATOR CIRCUIT SCHEMATIC 




Figure 1. 



Figure 2. Accelerator Assembly 




Figure 3. Current Sheet Properties 



Figure 4. Thruster on Test Stand 
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NONLINEAR TIME-DEPENDENT PHENOMENA 
IN REAL J X B DEVICES 

by 

S. T. Demetriades and G, S, Argyropoulos 
STD Research Corporation, Pasadena, California 



Even today - after 10 years of intensive effort — there is little agree- 
ment between theory and experiment in J X B devices. Previous theories 
do not provide a realistic description of the basic phenomena, because they 
make highly unrealistic oversimplifications and ignore physical mechanisms 
that have controlling influence on the performance and stability of these devices. 

Many interdependent physical mechanisms act simultaneously upon the 
flowing plasma and thus affect the local behavior (both in time and space) as 
well as the overall characteristics of a J X B device. These mechanisms 
and the resulting non-linear phenomena have become the object of realistic 
study at STD Research Corporation. This study has up to now analyzed the 
influence of thermal diffusion, thermal and velocity boundary layers, finite 
reaction rates, and electron energy convection on the local behavior and the 
overall characteristics in the stable mode of operation; the present research 
extends this analysis to the unstable (time -dependent) mode of operation. 

In MHD channels operating at high values of the Hall parameter cot, 
the fields and plasma properties undergo time -dependent fluctuations that 
are characterized by five different time scales. (By time scales we mean 
the characteristic period of such fluctuations. ) These time scales differ 
from each other by many orders of magnitude, depending on the geometry 
and the operating conditions. The overall performance of the device depends 
critically on which time scale is larger. If a physical process has a time 
scale that is orders of magnitude larger than the time scale of all others, 
temporal variations will be controlled exclusively by the slower process, 
while the faster processes can safely be assumed to proceed instantaneously. 

On the other hand, when such scales, arising from different physical 
mechanisms, are comparably large, then the performance depends on the 
interaction of all corresponding mechanisms. 

The five important time scales in MHD channels with segmented 
electrodes are : 



(1) The *4ime of flight** of the gas past an electrode period. 



tf = L/U, (1) 

where L is the electrode period and U is the gas velocity in the core of 
the flow; it reflects the influence of the geometry and of the boundaries on the 
time -development of the current and plasma property distributions, 

(2) The ionization -re combination time^ 



t 



R 



0. 433 / 





o 



(2) 



where r(Tg q) is the recombination rate for the three -body coll'isional 
recombination of the main ionizing species, evaluated at the reference tem- 
perature Tg Q, and n^ ^ is the electron number density corresponding to 
Saha equilibrium at Tg 
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( 3 ) 



(3) 



The electron energy- 



relaxation time 



1 



t 



E “ 





o 



T - T 
e, o g 





where Tg is the gas temperature, is the total collision frequency for 
electrons, is effective energy loss factor"^ for collisional energy 

transfer between electrons and heavy particles and y = 2kT where k 

is Boltzmann’s constant and* is the ionization potential of the main 
ionizing species. 

(4) The electron momentum transfer time^ 



t 

e 




(4) 



It is determined by the electron inertia term and describes the time necessary 
for acceleratioju of the electrons by the electric field to a given value of the 
drift velocity in the presence of collisions. 

(5) The time constant for electromagnetic propagation 



t = L/c , (5) 

em ' • ' ' 

where c is the velocity of light; it is defined from Maxwell’s equations. 

It is clear that t^^^ will be very short, and the same is true for t^, 
compared to the slower time scales t£, tp and t^*. Among the last three, 
it is not possible to distinguish one that will be the slowest under all operating 
conditions (see Table). All three, and the corresponding physical processes, 
must be considered simultaneously. This is done in our study, in contrast 
to previous investigations which'considered only t^,. 

This work has immediate application in evaluating the results of 
pulsed experiments which are used as models for continuous devices. Pulsed 
experiments are not subject to time scales that ^re bigger than the duration 
of the pulse, whereas continuous experiments are subject to all five time 
scales. Since each time scale is a function of the geometry and a very 
strong function of the operating conditions, their relative magnitude can be 
manipulated by the designer. Consequently, when the mechanisms that 
account for the difference between pulsed and continuous experiments have 
been identified, it is possible to define regimes by which continuous devices 
are governed by the same time scales as pulsed models. 

This work is based on a realistic numerical solution for the time- 
dependent current distribution. Such realistic solution is indispendable 
in a non-linear problem of this complexity, and it has shown that, under 
the influence of the aforementioned effects, the current distribution in 
J X B channels can indeed be very different from the predictions of 
previous theories. 

This work is complemented by experimental measurement of the time 
scale associated with the electron temperature rise in plasmas. Laboratory 
experiments consist of observation, by spectroscopic means, of the relaxa- 
tion of electron-temperature dependent properties of plasmas after a 
sudden change of electron temperature is produced by an externally applied 
microwave field, which heats the electrons but not (directly) the heavy 
particles, and therefore can selectively simulate some of the processes in 
J X B channels. This technique has also been used successfully to obtain 
laser output modulation through control of the electron temperature. 
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TABLE 

TIME SCALES IN TWO EXPERIMENTAL DEVICES 







Device 1 


- 




Device 2 




Gas: 


Nitrog 


en + 1 % Potassium 


Air 


+ 0. 2% Potassium 




L 




15. 0 mm 






42. 8 mm 




B 




2. 0 Wb/m^ 






2. 85 Wb/m^ 




P 




0. 5 atm 






4. 2 atm 




T = 

g 

U = 




3000 °K 






3054 °K 






3000 m/ s 






2925 m/ s 




t = 

em 

^f ~ 




1. 5 X 10‘^° sec 
5. 0 X 10"^ sec 






4. 3 X 10'^° sec 
1. 5 X 10"^ sec 




T 

e 


t 

e 


^R 


^E 


t 

e 


^R 


^E 


(°K) 


(sec) 


(sec) 


(sec) 


(sec) 


(sec) 


(sec) 


3200 


9. 9X10 


6. 3X10"^ 


1. 7X10"^ 


1. 6X 


10"^^ 1.6X10"^ 


4. 7X10' 


3400 


8. 3 


4. 1 


1. 0 


1. 4 


9. 9X10"^ 


2. 5 


3600 


7. 1 


3. 0 


6. 9X10'^ 


1. 3 


6. 8 


1. 7 


3800 


6. 2 


2. 3 


5. 2 


1. 2 


4. 9 


1. 3 


4000 


5. 5 


2. 0 


4. 1 


1. 1 


3. 8 


1. 0 


4200 


5. 0 


1. 9 


3. 3 


1. 0 


3. 2 


8. 4X10" 


4400 


4. 7 


1. 9 


2. 8 


9.4X10"^^ 2.8 


7. 1 


4600 


4. 6 


2. 0 


2.4 


8. 9 


2. 6 


6. 1 
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ENERGY TRANSFER MECHANISMS IN A PLASMA JET 

Adriano C. Ducati 

Plasmadyne, a division of Geotel, Inc. 

Santa Ana, California 

This activity has been primarily concerned with the study of interactions be- 
tween MPD thrusters and their environment. Strong evidence has recently 
been obtained that a substantial fraction of the arc energy in an MPD thruster 
is expended outside of the thruster, with a significant portion imparted to the 
ambient gas. Resulting interaction effects must be understood and evaluated 
so that the potential of this type of thruster operating in space can be deter- 
mined. It is difficult to make measurements that distinguish between thrust 
produced by acceleration of the propellant and by acceleration of the ambi- 
ent gas by direct action of the arc and associated magnetic fields. Two ap- 
proaches that are intended to complement each other are being used. First, 
the flow pattern in the vicinity of the thruster is being studied to see how 
much ambient gas is drawn into the active region near the nozzle exit. Pat- 
terns are determined by making a large number of surveys with calibrated 
orifice probes connected to sensitive gauges in a temperature-controlled en- 
closure. An earlier program using cold flow showed the entrained rate near 
the nozzle to be high. Measurements have now been extended to the hot gas 
case which, as would be predicted by theory, shows more rapid entrainment 
rates. Second, the behavior of the arc in the external region is being inves- 
tigated to establish mechanisms by which directed energy can be imparted 
to the ambient gas and to determine the importance of these effects. To ac- 
complish this, a thruster has been constructed for which the arc path is es- 
sentially all external. The complete arc region is accessible to probes and 
to photographic observation. Any thrust produced in this case is clearly a 
result of externally produced interaction forces. The unit can be operated 
with and without energizing a field coil, which permits results to be com- 
pared for cases with a magnetic nozzle effect as well as for cases for which 
self-field forces predominate. During the first tests, difficulties in starting 
and maintaining the stability of this unconventional geometry were encountered. 
A strong influence of the electrodes condition, especially of the cathode, was 
found. Another problem is concerned with the need for uniform erosion of 
electrodes and insulating parts so that a suitable geometry can be maintained 
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while these parts are continuously renewed by feeding. Techniques have re- 
cently been developed for taking motion pictures of the electrodes while the 
arc is in operation and with erosion taking place. These are very helpful in 
understanding the erosion process. They show molten material forming in 
drops that tend to concentrate the arc action and increase the local melting 
until the drop grows larger and is carried away by the flow. They also show 
the formation of cracks that tend to concentrate the arc, thereby increasing 
the thermal stress and the crack growth rate. Two newer areas of study are 
beginning to receive attention. First, the eventual interference of the opera- 
tion of a plasma jet with radio communication. A plume of ionized gas can 
refract or attenuate radio signals, and the operation of the arc can generate 
radio frequency noise that could mask weak incoming signals. Since the 
maintenance of communication is of serious concern, it is necessary to assess 
the importance of factors that may interfere with communication. The large 
fiberglass reinforced vacuum chamber in use at Plasmadyne is well-suited 
for a study of this type. Measurements that have already been reported sug- 
gest a correlation between the frequency of arc oscillations and the dominant 
frequency of noise produced. Work is now being initiated for determining 
the noise spectrum with different arc geometries and for evaluating possible 
ways of reducing the noise or of avoiding interference between the communi- 
cation frequency and the noise frequency. The second area of study is based 
upon the assumption that forces due to interactions with the ambient gas can 
be put to use to produce thrust in the upper atmosphere. Assuming the avail- 
ability of a sustained electrical power supply, a propulsion device using am- 
bient gas could produce a usable reaction force. A preliminary examination 
indicates that sufficient electrical power might be made available to over- 
come vehicle drag at altitudes above 200 kilometers. It is necessary to ex- 
plore the feasibility of influencing a large volume of rarefied gas in the free 
molecular flow regime, and to evaluate the possibility of sustaining an elec- 
trical discharge under these conditions. Although it seems difficult to simu- 
late the correct velocity and density conditions in the laboratory, useful 
information is expected to be obtained by testing MPD thrusters designed 
to produce large interactions with the ambient gas. Some of the tests made 
to determine laboratory measurement errors due to interaction effects can 
help to study the feasibility of the before- mentioned application. 

(Contract F44620-70-C-0002) 
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NON-EQUILIBRIUM PHENOMENA IN PLOWING PLASMAS 



M. Mitchner, C.H. Kruger, and R.H. Eustis 



In what follows, we describe two investigations currently 
in progress under this research program. References to recent 
publications of the results of other studies are listed at the 
end of this report. 



Non-equilibrium in Rapid Expansions 

Over the past decade, a considerable amount of study has 
been devoted to the phenomenon of vibrational relaxation of a 
single species of diatomic molecule. Characteristic times for 
the response of the vibrational temperature to either a sudden 
Increase (e.g., a shock wave) or a sudden decrease (e.g., a ra- 
pid expansion) in translational temperature have been measured 
and calculated. These processes are governed by the exchange of 
energy upon collisions, between the vibrational and the transla- 
tional modes, referred to as V-T exchanges. One of the outstan- 
ding problems in this area is that for some species the measured 
relaxation times for expansions are an order of magnitude smaller 
than for shock waves. 

In comparison with the preceedlng situation, very little 
work has been done on the relaxation of mixtures of different 
molecular species. In a non-equilibrium gas mixture containing 
more than one diatomic species, the vibrational temperature of 
each species can be different and the vibrational energies of 
the individual species can change through collisions Involving 
primarily an interchange in vibrational states, referred to as 
a V-V process. For polyatomic species, different vibrational 
temperatures need to be considered for each normal mode. 

A theory of vibration-vibration energy exchange between 
different species was first developed in 1952 by Schwartz, 
Slawsky, and Herzfeld [l]. However, no experimental data were 
reported until 1966. Taylor, Camac, and Peinberg [2] measured 
the V-V exchange between various molecular species behind a nor- 
mal shock wave. Measurements of V-V coupling for temperatures 
up to 1000°K have also been made using a laser-fluorescence tech- 
nique [ 3 ,^, 5 ] > For most mixtures, the Schwartz, Slawsky, and 
Herzfeld theory is in qualitative agreement with the experimental 
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data, but does not yield the correct magnitude nor the exact tem- 
perature dependence. Thus, the theory cannot be used to give 
quantitative results and experimental measurements are needed to 
obtain V-V coupling rates over a wide range of conditions. 

The purpose of this study is to investigate the vibrational 
Interactions of molecular mixtures using rapid expansions. In 
particular, the objective of the present effort is directed at 
investigating the V-V coupling between nitrogen and carbon dioxide. 
Data for this interaction will contribute to a quantitative un- 
derstanding of the operation of COg - Ng lasers. 

Preliminary measurements have been made of the nitrogen 
vibrational temperature in expansions of dilute N2 in argon 
and mixtures of N2 and CO2 in argon. These measurements were 
made using sodium D-llne reversal and spectroscopic line emis- 
sion methods. The measurements were made in a free jet. The 
use of the free jet enables measurements to be made as a function 
of distance in the flow direction without the necessity of con- 
structing separate viewing windows for each axial location. 

The results without CO2 are in agreement with those of 
Hurle, Russo, and Hall [6,7] in that they can be correlated 
using a relaxation time l/l 5 th of that measured behind a shock 
wave. Measurements for N2-CO2-A mixtures show that a small a- 
mount of CO2 has a large effect in decreasing the relaxation 
time of the nitrogen vibrational temperature. These data have 
been interpreted using a set of simplified rate equations, and 
indicate that the N2-CO2 V-V coupling is considerably stronger 
in an expansion flow than behind a normal shock wave. This pre- 
liminary result is the first indication that V-V relaxations may 
be faster in expansions than in shocks, in analogy with the situ- 
ation for V-T exchanges . 

Present work is directed at refining and supplementing our 
preliminary measurements. We plan to measure directly the relax- 
ation of the assymetric stretching mode of the CO2 using an infra- 
red band reversal technique. This measurement, in conjunction 
with measurements of the N2 vibrational temperature, will enable 
us to determine directly the N2-CO2 V-V interaction coefficient. 

The measurements will be compared with calculations based on a 
more general set of coupled rate equations. 

Ionization Instabilities 



The performance of non-equilibrium noble gas MHD generators 
has been found to be seriously affected by the presence of the 
ionization, or electrothermal, instability. The instability is 
observed as a fluctuation in the electron number density and all 
associated plasma properties. It has been found that there is a 
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critical value (Pcrit) Hall parameter cut at which the in- 

stability arises. Further increases in cot result in a c^ecrease in 
the effective conc3uctlvlty anc3 a saturation value Peff 
effective Hall parameter^ i,e. the arctan of the angle between the 
average current vector anc3 the average electric fielc3 vector in 
the gas frame of reference . 

It was suggestec3 by Brederlow and Hodgson [8] that Peff 
and Pcrit influenced by the size of the testing channel. 

Furthermore, Belousov [9] has indicated a dependence of wave- 
length on channel size. Gilpin [lO], on the other hand, predic- 
ted a wavelength independently of the channel boundaries. Ex- 
periments to more precisely determine the nature of channel- 
size effects are currently in progress. The material presented 
here describes the initial stages of this investigation. 

This portion of our program consists of an experimental 
study of the ionization instability in a continuous-electrode 
simulated MHD generator. The initial effort centered mainly on 
whether the behavior of the instability in a continuous-electrode 
experiment is consistent with the predictions of the present 
theories and with the existing experimental data taken in axial- 
discharge and Faraday-generator experiments. Our more general 
objective is to understand the effect of boundary conditions and 
channel size on the wavelength and amplitude of the instability 
and on the critical Hall parameter, the effective Hall angle, 
and the effective conductivity. Such information is expected to 
have an Important bearing on the extrapolation of results from 
present small-scale experiments to future MHD applications. 

Experiments have been conducted in atmospheric-pressure 
argon at 2300°K seeded with 0,11 % potassium. The channel had 
an electrode length of 15 cm and an electrode separation of 2cm. 

The gas velocity was 130 m/sec, the magnetic field varied up to 
1 Tesla, and the applied electric fields varied up to 500 v/m 
discounting the electrode losses. Measurements of the fluctua- 
tions in electric field were taken from a 1cm x 1cm region near 
the center of the electrodes by platinum probes projecting into 
the gas flow and separated by 5mm. The current distribution a- 
long each electrode was measured by segmenting the electrodes. 

It was found that the Instability followed closely the be- 
havior noted by previous experimenters and the predictions of 
the theories. At current densities less than 0.1 amp/cm^ and 
for Hall parameters cut up to 10 the instability was not apparent 
to any significant degree. At larger current densities the am- 
plitude and RMS value of the electric field fluctuations increased 
with current density and magnetic field. The critical Hall para- 
meter was found to be near 2. As cut was increased above 2 the 
effective Hall angle remained near arctan 2 and the fluctuating 
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component of the electric field Increased to 200 v/m (RMS). 

The regular fluctuations in electric field observed for cut = 4 
became increasingly Irregular at larger cut, suggesting breakup 
of plane wave structure. 

Voltage-probe data taken near the entrance region of the 
test section was much more irregular than that described above 
taken near the center of the section. Motion pictures taken in 
the entrance region show a very definite ionization front in- 
clined to the electrodes at an angle very near the effective 
Hall angle. 

In an attempt to short the fluctuating electric fields on 
a length scale much smaller than the electrode separation, three 
grids formed of l6 per inch tungsten wires were equally spaced 
between the electrodes. These grids were parallel to the elec- 
trodes and extended over the entire length of the electrodes. 

It was found that with the grids the instability was still pre- 
sent but the amplitude of the electric field fluctuations de- 
creased . 



More detailed experiments concerning the influence of chan- 
nel size on the wavelength of the instability and shorting effects 
in Faraday generators are currently in progress. 
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INVESTIGATION OF THE PRODUCTION OF 
HIGH DENSITY STABLE PLASMAS 

by 

R. M. Patrick, E. R. Pugh and D. H. Douglas-Hamilton 
Avco Everett Research Laboratory, Everett, Massachusetts 

I. INTRODUCTION 

A stable ionized plasma is required for many of the most interest- 
ing applications of plasma physics. 

Classically, stable plasma discharges are produced by working in 
the region where diffusion of the charged species to the walls stabilizes the 
discharge. The electron temperature in the discharge is raised (by adjust- 
ing the electric field E/P) until the ionization rate due to these electrons 
equals the rate at which ionized molecules are lost at the walls. However, 
the diffusion time becomes large compared to other instability times in the 
plasma, and the discharge spokes or arcs. Many of these instabilities are 
associated with the ionization process and its coupling with the electron tem- 
perature. AERL has constructed a facility under a company- sponsored pro- 
gram which produces uniform discharges in which the electron density, n^, 
can be varied independent of the electron temperature T^. Separation of the 
ionization process from the electron temperature is achieved in this facility 
by having a separate control on each, resulting in a stable discharge over a 
wide range of pressure and scale size. 

We achieve the separation of plasma ionization from the electron tem- 
perature by using a beam of high energy electrons to create a uniform ioni- 
zation throughout the volume of the plasma. Using the high energy electrons 
to provide the ionization, the electron temperature is varied independently 
by adjusting the electric field. It is found that if the electron temperature 
is kept low so that the ionization due to the discharge electrons is small com- 
pared to the ionization due to the high energy electrons, the discharge remains 
uniform over the range of pressures investigated, up to one atmosphere. As 
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the electric field is raised, the two ionization processes became equal and 
the discharge becomes nonuniform and eventually arcs. 

In the present program, we are using this facility to investigate the 
stability limits of such a discharge. These limits, of course, reflect upon 
the volumetric stability of other discharges important to practical applica- 
tions, such as those occurring in an extrathermal MHD generator. Volu- 
metric instabilities depend upon the recombination processes. This report 
describes some initial measurements at high pressure (300 Torr). Sub- 
sequent work will involve investigation of various discharge stability limits, 
relating these limits to specific applications. 
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II. DESCRIPTION OF AERL E-BEAM PLASMA FACILITY 



Theory 

Electron beam current is limited by the heating effect on the foil. 

The foil is not permitted to rise in temperature to a level where its struc- 
tural strength is significantly reduced, which would lead to pinholes and 
tears. The condition is that the total energy deposited in the foil during 
the pulse must produce a temperature rise less than ZOO^K. 

The fraction of the electron beam lost in the foil has been investi- 
gated by Borries and ‘Knoll, ^ and the rate of energy loss going through the 

2 

foil by Bethe and tabulated by Berger and Seltzer. For the aluminum foil, 
thickness 10“^ cm, and incident electron energy, 50 kV, used in the present 
experiments, the fraction of the incident beam energy lost in the foil is 80%. 



Assuming all the energy is deposited in the foil, a lower limit on the 
incident current density I (amp/cm ) in terms of the incident beam energy 
E (volts) and the pulse length (t) can be obtained. 

Elt < 0. 5 joule 

If E = 50 kV and t = 20 |j.sec, I < 0. 5 amp/cm . Approximately 20% of the 
incident electrons are transmitted, and the mean energy is reduced by 10 kV. 
These electrons are available for ionization of the gas. 

3 

The production rate of ions in a gas per cm is 



P = 



Ip < 
e 



8E 

9m 

E. 

1 



> 



-3 -1 

cm sec , 



where 

p = gas density, gm cm 

-2 

I = emergent electron beam current density, amp cm , 

9E 2 - 1 

<- Q^ > = mean energy loss rate, volt cm gm , 



^Borries and Knoll, Phys. Z. 35 , 279 (1934). 

^Berger and Seltzer, Tables of Energy Losses and Ranges of Electrons and 
Positrons, National Bureau of Standards. 
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= mean energy loss per ionization, volt, 
e = charge on electron, coulomb. 

The electron density n^ is given by 



dn 



dt 



e 2 

= -a n + p 
e ^ 



( 1 ) 



where a is effective recombination rate and p is production rate. Hence, 
in equilibrium, (times greater than the characteristic time T = l/an ), 
(dn^/ dt = 0), and 



n = 
e 



where 



aT 



-2 



Mm 

E. 

ke E. 

1 



( 2 ) 



P gas pressure, dynes cm 
T = gas temperature, ^K, 

3 -1 

a = effective recombination coefficient, cm sec , 

M = molecular weight of gas, 
m^ = proton mass, gm, 
k = Boltzmann's constant, erg/^K. 

Approximate maximum values of n for a typical electron beam, current 

2 ^ 

density 1 mA/cm , in N 2 and He are ^iven below, with characteristic decay 
time r and range R at E = 50 kV, using E^ = 50 volts. 
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When the electron beam is turned off, p = 0, and Eq. (1) may be 



integrated, giving 



n 

e 



n 



0 



1 + a n^t 



(3) 
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where Hq is the electron density at t =0. 

For values of time t » 1/ an^, we have 

Hg = l/at. (4) 

Thus, a value of a can be determined by measuring the value of n^ at 
time s t » 1/ a tlq . 

A pparatus 

A broad, fairly uniform beam of high energy electrons is generated 
by the apparatus shown in Fig. 1 and schematically in Fig. 2. Thermionic 
electrons from a tungsten filament are modulated by a grid whose potential 
can be varied with respect to the filament, and accelerated through a poten- 
tial Vq ( 50 kV) to a thin foil, in our experiments of aluminum. The foil 
divides the electron gun, maintained at a vacuum {< 0. 1 p), from the dis- 
charge chamber which may be at any pressure from below 1 p up to one 
atmosphere. 

After passing through the foil, the beam penetrates the discharge 
chamber through a wide mesh grid of stainless steel. This is used as the 
cathode, and a gold-plated disk opposite acts as the anode. Between these 
the plasma discharge voltage (< 10’,kV) is maintained. The grid is 
required to prevent damage to the foil, and the gold plate on the anode serves 

i! 

to reflect a proportion of the incident primary electrons, thus increasing the 
ionization in the gas. The grid structure containing the filaments is main- 
tained at -50 kV with respect to the foil (which is near ground) by a 5 pF 
capacitor which supplies the electron beam current. The filaments are 
pulsed negative 500 V with respect to the grid. It is estimated theoretically 
that approximately 1% of the electron beam current incident on the foil sup- 
port structure penetrates the foil and the cathode grid and is thus available 
for ionizing the gas. The foil support structure is a perforated plate to 
take the pressure of the gas in the discharge region. The effects of secon- 
dary electrons from these structures and ionization of the residual gas make 
this a difficult measurement, and it has not been done in a convincing manner 
as yet. 
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GAS SUPPLY 




The plasma discharge current is supplied by a 2 50 pF capacitor at 
voltages up to 10 kV. Either the anode or the cathode can be grounded. 

The gas being studied flows through the discharge chamber normal 
to the electron beam. Its velocity can be varied up to Mach 1, but up to 
the present time only quasistatic velocities (~ 1 M/s ) have been used, 
and only nitrogen has been examined. 

The uniformity of the electron beam in the discharge region and 
its approximate intensity variation was determined by replacing the anode 
wall with a lucite wall coated with sodium salicylate, a substance that is 
fluorescent when excited by high energy electrons. This was photographed 
and the intensity distribution shown in Fig. 3 was obtained. 



168 



I 




Q iq O • 

A1ISN3J.MI '/'ilVaa NObU.0313 03Znvt'!30N 



16 “ 



DISTANCE ALONG ANODE (cm) 



HI. INITIAL EXPERIMENTAL RESULTS 



Figure 4 shows a typical set of oscillogram traces showing the 
current and voltage measurements in the plasma discharge. The measure- 
ments were made at 300 Torr in N^ with 5000 volts across the plasma dis- 
charge with the anode grounded. The electron beam energy hitting the foil 
was 50 kV. The electron beam current between the pulsed grid and the 
foil is measured by an induction coil and is shown in the top trace of Fig. 4. 
The current in the discharge was measured across a resistor and is shown 
at two different gains in the next two traces. The potential distribution was 
determined by two probes 1. 5 mm in dia and 40 mm long placed in the dis- 
charge with a 10 Mfi impedance anode. Each probe was placed two cm 
from the anode or cathode and thus one cm apart. The voltages on the 
cathode and probes as well as the voltage difference between the two probes 
are also shown in Fig. 4. 

The current density was determined by isolating a region of the 

anode 2. 5 cm in dia and measuring the current to it. These measurements 

2 

indicated an effective anode area of 8Q cm . 



The time variation of E/P in a typical experiment (V^ = 5000 V, 

P = 300 Torr) is shown in Fig. 5. The corresponding value of V is 

.34 ^ 

obtained from Phelps * ; over the range of interest, it can be approxi- 

mated by: 



and 



Vj^ = 7. 5 X 10^ (E/P)2 
= f(E/P). 



The electron number density is determined from the current mea- 
surements and the drift velocity by the equation 

Engelhardt, Phelps and Risk, ‘'Determination of Momentum Transfer and 
Inelastic Collision Cross Sections for Electrons in Nitrogen using Transport 
^Coefficients, " Phys. -Rev. J_35, No. 6A, 1566(1964). 

Phelps, "Rotational and Vibrational Excitation of Molecules by Low Energy 
Electrons, " Westinghouse Research Laboratories, Scientific Paper 67-1E2- 
GASES-P2, 1967. 
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FIG. 4 Oscilloscope traces. 300 torr nitroge 
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FIG. 5. Electric field vs. time. 
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Ae V 



n 

e 



D 



where I is the discharge current shown in Fig. 4, A is the effective area 

of the discharge (80 cm^) and e is the electron charge. Values of n for 

6 

different times are shown in Fig. 6. The straight line has a slope -1 

indicating a l/t time dependence. The recombination coefficient given by 

a = thus 0. 8 X 10 cm^ sec* , in rough agreement with the 

results of previous workers for T ~ 0. 5^K. 

e 

11 - 3 

The maximum value of n^ obtained,*- 8X10 cm ,is in agreement 
with estimates made from the beam intensity and range in N^. 
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FIG. 6. Electron nvimber density vs. time. 
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IV. SUMMARY AND FUTURE PLANS 



AERL has a facility for making uniform discharges with variable 

at high pressures and large volumes. Presently, in this facility the 

discharge is uniform at pressures as high as one atmosphere with dis- 

2 

charge current densities up to 0. 1 amp/cm . 

Using this facility, some preliminary measurements were made 
concerning recombination rates at high gas pressure and variable T^. 

Initial measurements of the discharge voltage and current indicate that 

-7 -3 -1 

the effective recombination rate in N 2 is approximately 10 cm sec 
in agreement with previous work. 

Future work will validate these preliminary results and expand 
them to wider ranges of pressure, electron currents, discharge currents, 
and different gases. The limits on electron temperature for maintaining 

t 

a viniform discharge and the basic instabilities involved will be studied. 
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SUMMARY OF RESEARCH 



Supported by AFOSR Contract No. F44620- 69-C-0109 

J. Rosciszewski 

CALCULATIONS OF FLOW IN MAGNETIC NOZZLES AND FLOWS 
WITH RADIATION ENERGY EXCHANGE 

Continuation of research was conducted in three major groups: 

1. Problems associated with operation of an MPD arc-jet with 
pure magnetic nozzle. 

2. Transient flows with strong interaction with electromagnetic field 
or with intense radiation energy absorption. 

3. Explosion with radiation cooling wave. 

1 . MPD arc-jet with ”pure^’ magnetic nozzle . 

The first group pf problems is a continuation of research aimed toward 
understanding the acceleration process in an MPD arc-jet. This research was 
designed as a joint effort with A. Ducati' s experimental group at Plasmadyne 
Corp. In order to cut down the anode potential drop a device with a very 
short electrode was selected where the magnetic field is perpendicular to 
the anode and forms a ’’virtual electrode"(Fig. l)Because electrons can move 
easily along magnetic field lines, the anode drop should be reduced to zero 
and the device will operate as a pure magnetic nozzle. 

Experiments by A. Ducati indicated trouble by starting an arc in such 
a wall-less device. This is due to a fact that neutral gas exposed to a vacuum 
is streaming radially. The motion of neutral gas across an arc sheet, taking 
into account finite rate of ionization, was calculated on the computer and it was 
found that for /e^/^the current sheet thickness of this sheet should be at least 
j to prevent escape of neutrals. In order to cut down the escape of neutrals 
one could provide outside walls parallel to the magnetic nozzle (Fig. 1). 
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Figure 



Pure magnetic nozzle MPD arc-jet. 
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2. Transient flows with strong interaction with electromagnetic field or with 



intense radiation energy absorption. 

Last year's effort of simulation of non-steady shock-induced MHD flow 
is continued with alternation of a basic mathematical tool and different appli- 
cations. A modified Lax-Wendroff method was applied to shock-induced 
MHD flow in order to simulate experimental conditions of B. Zauderer. The 
resulting paper was presented last year at the Seventh International Shock ' 

Tube Symposium in Toronto, and will appear in transactions of the Symposium. 
General agreement with experiment turns out to be fairly good. However / there 
are some details which the present mathematical model cannot handle. One of 
the problems is difficulty in introduction of the relaxation time. This is due 
to the smeared shock thickness over a few mesh points. The present effort 

is aimed to application of the Lax-Wendroff method to equations in a system 

✓ 

of coordinates moving with the shock. In this way shock will be represented 
by discontinuity. At the same time one could obtain a secondary "smeared" 
shock resulting from strong interaction. This will make it possible to intro- 
duce relaxation time in the MHD problem. 

As an example of application of the method the flow behind a strong 
shock wave with absorption of a parallel beam of laser light is considered. 

3. Explosion with radiation cooling wave . 

The latest analytical effort is devoted to analytical treatment of 
transient flows with radiation energy exchange. The author's interest was 
induced by many still unsolved basic problems of high temperature flow 
important in the defense field. A classical problem is that of point explosion. 
Solution of this problem neglecting ambient pressure and assuming atmosphere 
with a uniform density was first given by Sedov (ref. 1) and by Taylor (ref. 2). 
This solution was obtained by introduction of a similarity parameter 

as an only independent variable. Here, small P denotes radius 
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and (tj shock wave radius. In this way the system of partial differential 

equations reduces to ordinary differential equations which could be solved in a 



closed form. 

Since the time of the classical Sedov-Taylor solution, many works have 
appeared including ambient pressure, atmospheric density variation, etc. 

An entire book by Korobeinikov has been devoted to the explosion problem 
(ref. 3 ). In this book most of the problems are solved by assuming small 
deviation from a similarity solution leading to linearized equation. 

Another very simplified solution assuming that all the mass is con- 
centrated at the shock wave was given by Chernyi (ref. 4). A recent paper 
by Laumbach and Probstein (ref. 5 ) also gives a similar simplified approach 
to the solution of a problem assuming concentrated mass at the shock. More 
general methods are based on finite difference treatment. However, such 
methods are extremely time-consuming and are always associated with 
numerical instability problems. Moreover, such methods do not take advan- 
tage of the known solution of a simplified case where similarity solution is 
applicable. It seems that any numerical method which does not include such 
a ’’skeleton solution,” which chould be easily obtained in the simplified case, 
could not be too successful. Recognizing these facts, the author of this 
proposal spent considerable time in obtaining a solution of an explosion 
problem using a very general method of integral relations. Such a solution 
is not restricted by an assumption of concentrated mass at the shock, etc. 

The method was checked in the simplified case where a similarity solution 
is applicable. A simple solution was obtained using asymptotic methods, and 
agreement with exact theory is excellent. 

As the next step, the author applied a multi-strip integral method 
to an explosion problem, including propagation of a cooling wave. Analysis 
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of singularity gives the initial condition for a cooling wave propagation. It 
was shown that the initial velocity of a cooling wave is exactly equal to relative 
gas velocity with respect to a shock wave. This is a kind of Chapmann- Jouget 
condition. A cooling wave must be therefore initially tangential to a shock 
wave (Fig. 2). 
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Radiation cooling wave 




Figure %. Explosion with radiation cooling. 
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